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ABSTRACT 

We explore the Magellanic Stream (MS) using a Gaussian decomposition of the H I velocity profiles 
in the Leiden-Argentine-Bonn (LAB) all-sky H I survey. This decomposition exposes the MS to be 
composed of two filaments distinct both spatially (as first pointed out by Putman et al.) and in 
velocity. Using the velocity coherence of the filaments, one can be traced back to its origin in what 
we identify as the SouthEast H I Overdensity (SEHO) of the Large Magellanic Cloud (LMC), which 
includes 30 Doradus. Parts of the Leading Arm (LA) can also be traced back to the SEHO in velocity 
and position. Therefore, at least one-half of the trailing Stream and most of the LA originates in the 
LMC, contrary to previous assertions that both the MS and the LA originate in the Small Magellanic 
Cloud (SMC) and/or in the Magellanic Bridge. The two MS filaments show strong periodic, undulating 
spatial and velocity patterns that we speculate are an imprint of the LMC rotation curve. If true, 
then the drift rate of the Stream gas away from the Magellanic Clouds is ~49 km s~^ and the age 
of the MS is ~1.74 Gyr. The Staveley-Smith et al. high-resolution H I data of the LMC show gas 
outflows from supergiant shells in the SEHO that seem to be creating the LA and LMC filament of 
the MS. Blowout of LMC gas is an effect not previously accounted for but one that probably plays an 
important role in creating the MS and LA. 

Subject headings: Galaxies: interactions - Galaxies: kinematics and dynamics - Galaxies: Local Group 
- Galaxy: halo - Intergalactic Medium - Magellanic Clouds - Radio Lines: general 



1. INTRODUCTION 

Under the prevailing concordance cold dark matter cos- 
mology, large structures like galaxies form through hier- 
archical accretion and merging of dark matter subhalos 
(e.g.. White & Rees 1978; Davis et al. 1985; Navarro, 
Frenk & White 1996, 1997; Moore et al. 1999). While 
much of the merging took place at early times, the pro- 
cess of accretion onto large spiral galaxies, such as our 
Milky Way (MW), continues at a reduced rate until late 
times (Bullock & Johnston 2005). Disruption and ac- 
cretion of small galaxies gives rise to gaseous and stel- 
lar tidal streams that continue to orbit the accreting 
galaxy as fossil relics of the cannibalistic activity. Many 
striking examples of disruption around our Milky Way 
have been discovered in recent years: the Pal 5 stream 
(Odenkirchen et al. 2001; Grillmair & Dionatos 2006), 
the Sagittarius stream (e.g., Ibata et al. 2001; Newberg 
et al. 2002; Majewski et al. 2003); the Monoceros stream 
(e.g., Yanny et al. 2003); the orphan stream (Belokurov 
et al. 2006; Grillm air 2006a); and the anticenter stream 
( |Grillmair||2006bD . There are likely numerous more such 
streams of stripped debris remaining to be discovered. 

But the most prominent and earliest discovered stream 
is the Magellanic Stream (MS), which stretches over 
100° across the southern sky behind the Large and Small 
Magellanic Clouds (LMC and SMC). From the mid-1960s 
onwards there were man y efforts to de tect high- velocity 
H I clouds, but Wanni er fc Wrixon| ( |1972) were the first to 
recognize the large extent of what was to become k nown 
as the Magellanic Stream. Mathewson et al. ( 1974 ) asso- 



ciated the Stream with the Magellanic Clouds and more 
fully demonstrated its ~100° span. A more thorough 
history an d general review of the MS is g iven bvJP^ 
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man et al. (2003, hereafter P03) and Briins et aL] ([2005 
hereafter B05). Whereas, except for the MS, all the 
above mentioned streams are stellar and attributed to 
tidal forces for their origin, the MS is still only recog- 
nized as a gaseous feature^ and is one such structure 
whose origin is still debated. 

Due to their proximity to each other and as the most 
massive of the MW satellites, the Magellanic Clouds 
(MCs) have long been considered to have influenced each 
other as well as the growth and evolution of their host 
galaxy. Thus, the MW-LMC-SMC system is regarded 
as an important laboratory with which to study the for- 
mation, evolution, and interaction of galaxies and their 
stellar populations. To this end, extensive mappings of 
the MCs have recently been conducted at a number of 
wavelengths. Surveys have been made in H I emission at 
21-cm (HIPASS: Barnes et al. 2001; Briins et al. 2005; 
Staveley-Smith et al. 1997, 2003; Stanimirovic et al. 1999; 
and Muller et al. 2003), in molecular spectral CO lines 
(NANTEN: Fukui et al. 1999; Mizuno et al. 2001), as weh 
as in the radio continuum (Haynes et al. 1991; Dickel 
et al. 2005), in the thermal infrared (IRAS: Beichman 
et al. 1988; MSX: Mih et al. 1994; SAGE: Meixner 
et al. 2006; S^MC: Bolatto et al. 2007), in the near- 
infrared (DENIS: Epchtein et al. 1997; 2MASS: Skrutskie 
et al. 2006), in broadband optical colors (MCPS: Zarit- 
sky et al. 2002, 2004) as well as in optical emission lines 
(MCELS: Smith et al. 1998), and at ultraviolet (Smith, 
Cornett & Hill 1987) and X-ray wavelengths (ROSAT: 



^ There is a sparse stellar population at the location of the ] 
ellanic Bridge feature (Irwin, Demers &: Kunkel 1990), but their 
connection to the H I Bridge is unclear, as is the connection of the 
Bridge to the classical MS. 
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Snowden & Petre 1994; Chandra and ACIS: Townsley et 
al. 2006). While these surveys have deepened our un- 
derstanding of the Clouds themselves — e.g., their star 
formation history, their stellar content, and their overall 
structure (van der Marel 2001) — the H I investigations 
are central to understanding the most obvious product 
of the MW-MC interaction — the Magellanic Stream. 

Obviously, the overall appearance of the MS is shaped 
by the dynamics of the MCs wh ich are moving almost en- 
tirely tangentially in the sky (Kallivayalil et al. '2006ajbj 
[Piatek, Pryor, fc Olszewsk i 2007) and have recently just 
passe d perigalacticon in their polar orbit around the 
MW dBesla et aL]|2007D . The earliest dynamical stud- 
ies of the MCs used the MS to constrain the orbit of the 
LMC and thereby obtained an LMC s pace velocity of 
^350 km s~^ (M urai & Fujimoto 1980) |Lin fc Lynden- 
|Ben||1982| |HelleF& Rohlfs 1994- kowever, most proper 

CMC 



However, more recent large scale H I surveys have re- 
vealed new complexities in the MS that are difficult to 



motion measurements of the 



(summarized in 



van 



der Marel et al. |2QQ2, hereafter vdM02) have f avored 



a lower spa ce velocity of around 250 km s ^. [Murai 



& Fujimoto ( [l980} , and similar subsequent MS model- 



mg papers, used these space velocities and an isothermal 
sphere MW potential to derive orbits for the MCs with 
an orbital period ~1.5 Gyr and a last apogalacticon dis- 
tance of ~120 kpc. All of these models agree that the 
MCs had a close encounter ~200 Myr ago. However, the 



new HST proper motions of the MCs (Kal livayalil et al. 
2006a|b[ ) give MC space velocities ~100 km s ^ higher 
than those produced by the earlier proper motion sur- 
veys and that increase the orbital period to ~2-3 Gyr 
and increase the orbital elliptici ty (i.e. the last apogalac- 
ticon distance at ^ 150-200 kpc; [Kallivayalil et al.|2006b] 
Besla et al. 2007). But the addition of a more realistic 



NFW potential for the MW produces the startling result 
of hyperbolic orbits for the MCs (Besla et al. '20071) Even 
with a high MW mass model, while the MCs "become 
bound" again, the LMC orbital period is ~7 Gyr and the 
last apogalacticon distance is ~400 kpc. Clearly, whether 
the MCs are bound or not, and the shape of their orbits, 
has a direct influence on the interactions that produce 
and shape the MS. 

Large- area 21-cm radio surveys have produced most of 
the information now available about the detailed struc- 
ture of the MS. Since its discovery as a long stream of 
H I gas trailing the MCs a number of models have at- 
tempted to explain the dynamics and origin of the MS. 
Early N-body simulations with hundreds of particles by 
Lin & Lynden-Beh (1977, 1982) and by Murai & Fuji- 
moto (1980) were able to reproduce the general features 
of the Stream (such as its length and velocity distribu- 
tion) through tidal stripping by the MW. Later on, it 
was proposed that the MS could have been created by 
ram pressure forces (Meurer, Bicknell, & Gingold 1985; 
Moore & Davis 1994) as the MCs move through the hot 
gaseous halo of the MW. Ram pressure strips some gas 
from the Clouds and creates a trailing gaseous stream. A 
persistent problem with the tidal models is that they pre- 
dict a stellar MS component that, to date, has not been 
observed despite numerous efforts (e.g., Philip 1976a,b; 
Recillas-Cruz 1982; Briick & Hawkins 1983; Kunkel et 
al. 1997; Guhathakurta & Reitzel 1998). Because ram 
pressure only affects gas and not stars, these models 
seemed initially to be more consistent with a gas-only 
structure. 



account for in a r am pressure model. For example. Put 



man et al. ( [1998[ ) used the HIPASS data ([Barnes et aL 
200 ip to mscover a gaseous leading arm of the MS. This 



IS a feature readily accounted for by the tidal models 
but creating a leading arm by ram pressure forces re- 
mains a formidable problem. The HIPASS data have 
also shown that the trailing MS is spatially bifurcated 
(P03), although Cohen (1982) and Morras (1983) previ- 
ously pointed out that the MS splits into two branches. 



Mastropietro et al. (2005) performed a large ram pres- 
sure-ftidal force simulation of the Magellanic Stream 
(with only the LMC as a progenitor) and were able to 
reproduce the general features of the Stream, including 
its extent, shape, column density gradient, and veloc- 
ity gradient; however, Mastropietro et al. could not re- 
produce the spatial bifurcation of the Stream nor the 
Leading Arm Feature. The recent N-body tidal simu- 
lations by Connors et al. (2004, 2006), in which most 
particles are stripped from the SMC during a close en- 
counter with the LMC and MW ~1.5 Gyr ago, give the 
closest reproduction of the Stream to date, including the 
spatial bifurcation of the Stream, the Leading Arm (and 
its bent shape), and the MS velocity distribution. How- 
ever, a problem with most tidal models, including those 
by Connors et al., is that they have trouble reproducing 
the column density gradient along the Stream, whereas 
ram pressure models match this particular feature of the 
observations better, and, of course, account for an en- 
tirely gaseous Stream. While there is some evidence that 
stellar tails may be offset (or completely missing) from 
the gaseous tails in tidally interacting galaxies (Mihos 
2001; Hibbard, Vacca, & Yun 2000), the study of [John- 



ston 
shouL 



S1998) indicates that stellar debris from the LMC 
r have already been found. 
Most of the current literature has supported either the 
SMC or Magellanic Bridge as the source of the MS gas 
(e.g., P03, B05) since it appears to emanate from these 
regions in maps of H I column density on the sky. More- 
over, the mass of the SMC is much less than that of the 
LMC (by an order of magnitude) and it is therefore, pre- 
sumably, much easier for SMC gas to be stripped than 
LMC gas. Many of the tidal models have used the SMC- 
origin assumption in their N-body simulations (Gardiner 
& Noguchi 1996; Yoshizawa & Noguchi 2003; Connors 
et al. 2004, 2006). On the other hand, Mastropietro et 
al. were able to reproduce the general characteristics of 
the MS by ram pressure stripping from the LMC alone. 

In this paper we follow the tradition of using large- 
area 21-cm data to investigate the relation of the MS 
and the MCs. We take advantage of the high velocity 
resolut ion of the Leiden- Arg entine-Bonn (LAB) H I dat- 
acube (Kalberla et al. 2005) to investigate the detailed 
structure of both the leading and trailing arms of the MS 
across their entire known length. Particular attention is 
paid to disentangling MS features from other overlap- 
ping structures in the datacube, including the MW disk, 
other Intermediate- and High- Velocity Clouds, and the 
MCs themselves. This has allowed us to uncover several 
key aspects of the MS that differ from the earlier inter- 
pretations and models of the MS and that lead us to a 
new mechanism to explain the origin of the Stream. 
First, we show that one of the two trailing MS fila- 
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merits as well as the Leading Arm originate in the LMC, 
not in the SMC or in the Bridge as previously suggested. 
Moreover, the specific site that is the source of this gas 
we identify as the SouthEast H I Overdensity (SEHO), 
a region of dense H I and intense star-formation in the 
southeast of the LMC. Analysis of th e high spatial res- 
olu tion H I Parkes data of the LMC ( St aveley- Smith et 
al. 2003) indicates that supergiant shells in the SEHO 
are probably responsible for blowing out much of the gas 
from the LMC, creating (at least one of) the filaments 
of the MS and Leading Arm. Once blown out and free 
of the gravitational grip of the LMC, the gas experiences 
tidal stretching from the MW potential and separates 
into the leading and trailing components. This blowout 
mechanism for releasing Magellanic gas represents an al- 
ternative to the tidal and ram pressure models. Finally, 
we propose that the periodic sinusoidal weaving of the 
trailing Stream filament may be a result of the off-center 
position of the SEHO in the rotating LMC disk. Cou- 
pled with knowledge of the rotation period of the LMC 
at the SEHO radius, we can use the sinusoidal patterns 
to estimate the total age of the MS as ~1.74 Gyr under 
the proposed scenario. 

This paper is organized as follows: Section 2 is a de- 
scription of the 21-cm LAB data. In Sections 3 and 4 we 
describe the Gaussian Decomposition of the LAB data 
and the separation of the MS from the MW features; the 
casual reader uninterested in these details may wish to 
proceed to Section 5. In Section 5 we present the results 
of our investigation of the structure of the Magellanic 
Stream using our database of Gaussian Decomposition 
centers. Section 6 details our analysis of the SE H I 
overdensity in the LMC. A discussion of our findings and 
conclusions are given in Section 7, and a summary of the 
primary conclusions is given in Section 8. 

2. BRIEF DESCRIPTION OF LEIDEN-ARGENTINE-BONN 
(LAB) DATA 

The Leiden-Argent ine-Bonn all-sky H I survey 
Kalberla et al. 2005) is a combination of the Lei- 
en/Dwingeloo Survey (LDS: Hartmann & Burton 
1997), covering the sky north of S = —30°, and the Insti- 
tuto Argentino de Radioastronomia Survey (lAR: Arnal 
et al. 2000; Bajaja et al. 2005) at more southern declina- 
tions. The combined material has a velocity resolution of 
1.3 km s~^, a spatial resolution of 36' on a grid spacing 
of 0?5 in Galactic latitude {b) and 0?5/cos(6) in Galactic 
longitude (/), and is corrected for stray radiation. The 
velocity range of —450 km s~^ to +400 km s~^ is ade- 
quate for (almost) all Galactic work. The root mean 
square noise is 0.09 K. Here, we exploit the extensive sky 
coverage of the LAB data and its relatively high velocity 
resolution to follow the MS in detail over its full length, 
and to unravel its filaments. 

3. DESCRIPTION OF AUTOMATED GAUSSIAN 
DECOMPOSITION 

To improve our ability to trace structures of the MS 
and disentangle them from MW gas we wrote an auto- 
mated Gaussian analysis program in the Interactive Data 
Language (IDL)^ using an algorithm similar to that used 

^ A product of ITT Visual Information Systems, formerly Re- 
search Systems, Inc. 



by |Haud| ( |2QQQ^ ). Kalberla & Hand (2006) have also per- 
formed a Gaussian decomposition of the LAB database 
but with different goals, namely the physical interpreta- 
tion of the structure of High Velocity Clouds. Our use of 
Gaussian decomposition is predicated upon the expected 
continuity of the filamentary structures in terms of veloc- 
ity, position, velocity-dispersion, and integrated column 
density, which allows us to track features even through 
complex, crowded H I environments. 

The general algorithm to decompose an H I velocity 
profile into Gaussians proceeds in two stages. In the 
first stage new Gaussians are added to the model fit un- 
til the root mean square of the residuals (observed — 
fitted profile; rms) drops below the noise level, or any 
new Gaussians do not improve the model fit significantly 
{5rms <2%). In the second stage an attempt is made 
to reduce the number of Gaussians in the fit without in- 
creasing rms appreciably. The details of these stages are 
discussed in the following sections. 

One of the benefits of the Gaussian analysis performed 
here is that faint H I structures can be enhanced by plot- 
ting the integrated column density of each Gaussian at 
its central velocity. This avoids spreading the flux of the 
Gaussian over a range of velocities and instead concen- 
trates it to one point. Figures [4) (e) (g) [lO) [llj [M] and [17 
in this paper use this technique. 

3.1. Gaussian Fitting 

In order to find the best-fitting Gaussian decomposi- 
tion to a velocity distribution along a given line-of-sight 
we use the general purpose IDL curve-fitting package 
MPFIT written by Craig Markwardt.^ MPFIT is a set 
of routines for robust least-squares minimization (curve 
fitting) based on the MINPACK-1 FORTRAN package. 

In order to adopt MPFIT to our purposes we wrote 
an IDL function that returns the output of the Gaussian 
function, Tb = Tb,o exp( — (i; — vo)'^/2al)^ given an ar- 
ray of velocities and the Gaussian parameters (Tb, vq, 
and ay). The necessary inputs for MPFIT are the set of 
values for the independent and dependent variables (the 
observed data: velocity - and brightness temperature 
- Tb) and a first guess for the Gaussian parameters. One 
can also give upper and lower limits for each parameter 
that is allowed to vary in the fit; this is a useful feature 
for constraining the Gaussians to meaningful solutions. 
MPFIT returns the parameters for the best fit, the for- 
mal 1-cr uncertainties for those parameters, and of the 
fit. We also used rms (root mean square of the residuals) 
as another means to ascertain the "goodness" of the fit. 

Our Gaussian decomposition program uses MPFIT to 
fit Gaussians to peaks in the H I profiles or to peaks in 
the residuals (profile — previous best fit). For the ini- 
tial guess of the Gaussian parameters for a given single 
peak we use the height of the peak for Tb^o, and the 
velocity of the peak for vq. To get an initial guess for 
ay we used the fact that the ratio of the derivative of 
the Gaussian to itself is a line with a slope of — 1/cr^; 
Tb/Tb = (-1/cr^) X {v -vo). A line was fit to this ratio 
near the central part of the peak, and ay computed from 
its slope. In MPFIT the parameters were constrained 
to lie between the limits: 0.01 K < Tb,o < 2 x TB,max, 

^min < ^0 < ^max, 0.5 km S"^ < (J^ < ^ X (v^ax " ^min); 



Available at http://cow.physics.wisc.edu/~craigm/idl/idl.html 
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where Te^max is the maximum Tb of the whole profile, 
and '^min and ^max are the minimum and maximum ve- 
locities of the whole profile (—450 and +400 km s~^). 
Tb,o was also constrained so that only peaks higher than 
a certain threshhold a bove the noise level were chosen for 
Gaussian fitting (see ^3. 2| below). 

To obtain the noise level we smoothe d the H I profile 
with a [1 6,16,2] Savitzky-Golay filter (Savitzky & Go 
lay 1964|) (where the numbers in brackets refer respec- 
tively to the number of data points to the left and right 
of each point to include in the filter, and the order of 
the derivative desired) and then subtracted this from the 
original H I profile to remove any real features. Then we 
found the standard deviation of points with \v\ > 250 
km s~^ after 5a outliers were rejected. This was used as 
the noise level for the given H I profile. 

3.2. Adding New Gaussians 

Gaussians were added to the velocity profile fit one at 
a time. The current best fit was subtracted from the 
observed H I profile to find the residuals. These residu- 
als were then searched for peaks higher than 5 times the 
noise level (although this was lowered to 2 times the noise 
level if no peaks were found) and MPFIT was used to find 
the best-fitting Gaussian for each peak along with its rms 
using the first guesses described in the previous section. 
To smooth over the noise and search for features on var- 
ious scales five smoothed versions of the residuals were 
created using Savitzky-Golay filters of [4,4,2], [16,16,2], 
[30,30,2], [50,50,2] and [100,100,2] and Gaussians were fit 
to all the peaks in these smoothed profiles. 

The best-fit Gaussians to peaks in the six versions of 
the residuals (the original and the five smoothed) that 
had parameters within our acceptble limits (see above) 
were kept for further fitting. Each one of these candi- 
date Gaussians is taken in turn and added to the current 
best-fit Gaussian decomposition of the entire H I profile. 
MPFIT is then rerun with this new Gaussian decom- 
position as a first guess, fitting the multiple Gaussians 
at once, and the change of rms compared to the previ- 
ous best-fit was computed. The candidate Gaussian that 
gives the greatest decrease in rms is then added to the 
overall decomposition of the H I profile. 

This procedure is repeated and Gaussians added to the 
decomposition until the rms drops to or below the noise 
level, or the decrease in rms is less than 2%. 

3.3. Removing Gaussians 

The best-fitting Gaussian decomposition is that which 
minimizes both the rms as well as the number of Gaus- 
sians. To achieve this goal, we attempt to remove Gaus- 
sians that do not significantly improve the fits to the 
velocity profiles. The Gaussians in the best-fit decom- 
position are sorted in order of area {A = Tb,o x cTvV^ 
K km s~^), and the smallest half of the Gaussians were 
picked for possible removal. Each one in turn was tem- 
porarily removed from the decomposition and the new 
best-fit and rms for the whole profile are found with 
MPFIT. If the increase in rms is less than 2% the Gaus- 
sian is permanently removed from the decomposition. 
Even though it might appear at first glance that this 
is repeating work done in the "Adding Gaussians" stage, 
this particular decomposition might not have been looked 
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Fig. 1. — Examples of Gaussian decomposition for LAB H I 
profiles (top plot in each panel) and their residuals (bottom plot 
in each panel). (a) The Gaussian decomposition at (l,b) = 
(51.0°,— 80.0°), with four Gaussian components (each shown by a 
different color, their sum by red), showing the Magellanic Stream 
at negative velocities, (b) The Gaussian decomposition at (l,b) = 
(296.0°,— 58.5°), with five Gaussian components, showing two sep- 
arated Magellanic Stream components at positive velocities. 

at before because of the order in which Gaussians were 
added. Gaussians are also removed from the decompo- 
sition if two quite similar Gaussians were found at the 
same velocity. In that case they are replaced by a single 
Gauss i an wi th parameters given by Equations 11-13 in 
Hmidl (12000). 



3.4. Selecting the Next Position 

We initially used the same procedure as 'Hand ( 2000 ) 
(see his section 3.3) to select the next position on the sky 
to decompose. If the profile at a neighboring position has 
already been decomposed but has a worse decomposition 
(either larger rms or more Gaussians) than the decom- 
position at the current position, then the program tries 
to re-decompose the neighboring profile using the best- 
fit solution of the current position as the initial guess. 
The scheme also allows the program to wander around 
re-decomposing H I profiles until the re-decomposition is 
not an improvement. We found that this "wandering" 
scheme was too CPU intensive and did not improve the 
solutions substantially. Therefore, we used a modified 
scheme that forced the program to return to the previ- 
ous position after re-decomposing a neighboring position. 
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Fig. 2. — Average number of Gaussians in the H I decomposition 
per position as a function of Galactic latitude b. 

For new positions, that hadn't been decomposed yet, the 
best-fit solution at the previous position was always used 
for the initial guess. 

For low latitudes the program was not allowed to re- 
decompose profiles since it took much longer in these 
regions. Our program did not take into account any self- 
absorption, so some low latitude profiles near the Galac- 
tic center will not be correctly represented by the Gaus- 
sian decomposition. This does not affect our study of 
the MS here because it should not have significant self- 
absorption, if any. 

3.5. Statistics of the Gaussians 

Our automated Gaussian decomposition program was 
run on ah 259,920 H I profiles (720 values of / x 361 val- 
ues of 6, in steps of 0?5) of the LAB all-sky survey. In the 
end, the entire sky was decomposed into 1,370,801 Gaus- 
sians. Several examples of Gaussian decompositions at 
various sky positions are shown in Figure [l] We find 
that the majority of H I profiles (at latitudes above the 
disk, \b\ > 15°) are well fit with four to six Gaussians. 
The average number of Gaussians per profile as a func- 
tion of b is shown in Figure |2j the number peaks at 12 at 
the Galactic center and levels off to four near the poles. 
The majority of Gaussians at higher latitudes are from 
local MW, zero-velocity gas. The distribution of fitted 
Gaussian parameters for various populations are shown 
in Figure ID a "2D histo gram" indicating the number of 
Gaussians with a particular height (Tb,o) and Gaussian 
width (ay). An intriguing structure is apparent in the 
distribution of zero- velocity gas Gaussians in Figure [3}d 
following di trend (for 8 ^ > 35 km s~^ and 1 

^ ^B,o ^ 3 K) and nearly conserving its area at ~35 K 
km s~^. It is not clear what this structure corresponds 
to. The distributions of all Gaussians as a function of 6, 
and vi^sR are shown in Figure [4] The disk of the MW, the 
local zero-velocity MW gas, and the Magellanic Clouds 
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Fig. 3. — The distribution of Gaussian parameters (Tb and ay) 
for various H I populations, (a) All Gaussians; (b) local zero- 
velocity Gaussians (|Vlsr| < 30 km s"-*^); (c) Milky Way disk 
Gaussians with |VlsrI > 30 km s"-*^ (see §4.2 for how these are 
defined); (d) Magellanic Clouds and Stream Gaussians. Note the 
different color scalings for a+b and c+d. 

and Stream are readily apparent. 

3.6. Validity of the Gaussian Decomposition 

It might be asked whether Gaussian decomposition is 
the correct way to analyze these data. Though Gaus- 
sian decomposition has been widely used as a tool to 
analyze all forms of H I data going back to the 1960s 
(e.g., Kaper et al. 1966; Takakubo & van Woerden 1966; 
Burton 1970; Schwarz & van Woerden 1974; see review 
by Hand 2000), it is only physically well-motivated for 
isolated, internally virialized clouds, and here it finds its 
most common application (e.g., Briins, Kerp & Pagels 
2001; Wakker, Oosterloo, & Putman 2002; Kalberla & 
Hand 2006). In contrast, here we are using Gaussian 
decomposition primarily as a tool to disentangling over- 
lapping H I structures under the presumption that H I 
structures along the line-of-sight only slowly vary those 
properties encapsulated by a Gaussian description (ve- 
locity, position, velocity-dispersion, and integrated col- 
umn density). 

We have looked at various isolated H I clouds in the 
LAB data and found that they are well-fit by Gaussians. 
Moreover, by using Gaussians we are able to disentangle 
different H I filaments even when they are overlapping 
in velocity. In those situations it is clear that the Gaus- 
sian decomposition traces structures that are real and 
they may even hold physical information about the struc- 
tures. We are still successful in tracking tenuous struc- 
tures through even more complicated environments even 
though the decompositions of those environments (local 
MW zero- velocity gas, the MW disk, and the Magellanic 
Clouds) likely holds no physical meaning. Whether or 
not it holds physical meaning, we are primarily interested 
in how the Gaussian decomposition provides a represen- 
tation of the data that enables us to track large features 
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Fig. 4. — Position-velocity distribution of detected H I Gaussians in the LAB datacube where each is represented by a single point with 
weight equal to the Gaussian area and the results summed over all Gaussians. (a) Vlsr vs. I, summed along b, showing the characteristic 
velocity curve of the Milky Way as well as the Magellanic Clouds and Stream at larger velocities; (b) Vlsr vs. 6, summed along I, showing 
the zero- velocity gas (even at high latitudes), Milky Way disk (I ~0), Magellanic Clouds (to upper left) and Stream (arcing on left side). 
The greyscales are in units of K km s"-*^. 



of the MS. In addition the decomposition allows us to 
reduce our data to a manageable size (from a datacube 
to a database of Gaussians). 

There are a high number of Gaussians with low Tb,o 
and ay which are most likely due to fitting of noise fea- 
tures. We make a simple parameter space boundary 
^B,o > —0.08(7^ + 0.45 to remove these "noise" Gaus- 
sians. The presence of possibly extraneous narrow or 
wide Gaussians does not influence our interpretation of 
the LAB data, which is based on significant features in 
the database. 

4. REMOVAL OF MILKY WAY GAS 

In order to study the MS, we need to separate the 
Gaussians of MS gas from those of MW gas, which is a 
particulary difficult problem at low and at any lati- 
tude when the velocity of the MS differs little from the 
velocities expected for gas in the conventional MW disk. 
The velocities of the conventional MW disk gas roughly 
follow the expectations for material moving in circular 
orbits around the Galaxy, i.e. 



V^SR{l,b,R) 



Rq 

~R 



V{R) - Vo 



sin(0cos(6) (1) 



where Vq and Rq are the solar velocity and Galactocentric 
distance, respectively, and V{R) is the rotation curve at 
R. At high latitudes the gas is (mostly) concentrated to 
velocities of Vlsr ^ (zero- velocity gas) because the gas 
is (mostly) nearby, since we are looking out of the plane of 
the disk, and therefore Rq/R Essentially, this gas is 
local ISM gas moving in nearly the same orbit as the Sun 



around the MW. At lower latitudes the MW disk gas has 
a larger range of velocities of approximately, —Vocos{b) < 
"^LSR ^ +l^cos(6), but of course there is a strong sin(/) 
dependence. Due to the different characteristics of high- 
latitude versus low-latitude MW gas different methods 
were employed to separate them from the MS gas. It 
is worth emphasizing that ultimately no conclusions in 
this paper depend on the details or ultimate accuracy 
of these population decompositions; our care in pursuing 
these strategies is to make improved maps of the MS. 

4.1. Zero- Velocity Gas at High Latitudes 

At most positions in the sky, the zero- velocity MW disk 
gas is easily distinguishable from the MS gas because 
they have very different velocities. However, since the 
MS stretches from Flsr - 300 km s"^ to Flsr - -400 
km s~^ it must cross Vlsr = at some point. This hap- 
pens in the region -84° < b < -78° and 288° < / < 
327°. In this area it becomes challenging to distinguish 
the zero- velocity MW gas from the MS gas. It is difficult 
to disentangle the two populations of gas without some 
decomposition scheme as we have used here. Most earlier 
column density maps of the MS show gaps in this region 
(see Figs. 4-5 in P03, and Fig. 2 in BOS) which we intend 
to remedy in our maps. 

We attempted to separate the zero- velocity MW Gaus- 
sians from the low-velocity MS Gaussians by using the 
Gaussian parameters ay and Tb,o alone, since velocity 
would not be of much use in this case. However the zero- 
velocity MW and MS Gaussians also overlap in ay — Tb,o 
space (see Fig. [3| which makes this separation strategy 
untenable. The zero- velocity gas was eventually removed 
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Fig. 5. — A snapshot of the analysis used to remove the zero- 
velocity Milky Way Gaussians in the region where they overlap 
the Magellanic Stream. At a given Galactic latitude the central 
velocity of each Gaussian is plotted against its Galactic longitude. 
Each Gaussian is represented as a line, where its length corresonds 
to 2av and the color corresponds to Tb (the colorbar shows the 
scale). The continuity of features in position, velocity, Tb, and 
(Jv is used to identify them as due either to local Milky- Way zero- 
velocity gas or to the Magellanic Stream. An interactive program 
was developed whereby any Gaussian can be removed from the dat- 
acube by clicking on it on the display screen. The labels "STOP" , 
"ZOOM IN" , etc. are buttons on the interactive display to perform 
actions during the Gaussian de-selection process. 

interactively. We made Flsr vs. / plots of all of the Gaus- 
sians at a single h in the region where the MS crosses 
"^LSR = 0- Each Gaussian was represented by a verti- 
cal line, where the length of the line corresponded to the 
Gaussian's 2cr^ and the color of the line to its Tb,o- Fig- 
ure [5] shows an example of one of these plots, and demon- 
strates the relative ease with which such a representation 
makes it possible to distinguish the MS Gaussians from 
MW Gaussians because of the nearly-constant, but dif- 
ferent, Tb,o and trends for each and the straight ver- 
sus arcing trends that differentiate them. An interactive 
program allows us to remove any Gaussian represented 
in this way by clicking on it. All of the Gaussians con- 
sistent with being zero-velocity MW Gaussians were re- 
moved at a given 6, and the process was repeated for all 
h where the MS and zero- velocity MW gas overlapped. 
This process was iterated a few times to ensure that no 
residual zero- velocity patterns were left over. The re- 
sults of this zero-velocity MW Gaussian removal scheme 
can be seen in Figure |6] For the regions where there is 
no overlap between MS gas and Milky Way zero- velocity 
gas all Gaussians with |Vlsr| < 45 km s~^ have been 
removed. A separate scheme was used to remove the 
MW gas at somewhat higher velocities (i.e. |Vlsr| > 45 
km s~^ which is described in the next section. 

4.2. Milky Way Disk Gas at Low Latitudes 

Since we were interested in investigating the Leading 
Arm, we also needed to remove the MW gas at low lati- 
tudes, because the Leading Arm passes through 6 = 0°. 
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Fig. 6. — Vlsr vs. 6 of H I Gaussians (summed along I) (a) 
before the Milky Way and zero- velocity H I gas has been removed, 
and (b) after it has been removed. The greyscale is in units of K 
km s"-*^. 

We first attempted to use a simple, symmetric analytical 
model of the MW to remove the disk gas, but this failed 
to remove a significant portion of MW gas because of the 
exaggerated simplicity of the model. We instead adopted 
an empirical method to remove the disk gas. At a given 
/ and 1^SR7 the Tb due to the MW disk gas drops off 
quickly with b. The basis for our empirical strategy is 
to find where this drop-off occurs and call that the end 
of the disk. We used the profiles from the original H I 
datacube to accomplish this. In order to see the global 
trend and remove noise, each profile was smoothed with a 
Savitzky-Golay [15,15,2] filter. For each / and Vlsr > 30 
km s~^ the "edges" of the disk (where it fell below 5% of 
the central Tb value) in b were found. The values of these 
"edges" were then used to remove all Gaussians that fell 
within this b range at that particular / and Vlsr- In ad- 
dition, all Gaussians with Vlsr < 30 km s~^ have been 
automatically removed. The Gaussians left-over after the 
MW disk-gas removal (primarily HVCs, some IVCs, and 
MS gas) can be seen in Figure [7| The distributions of 
the Gaussian parameters (Tb,0 7 <^v) of the MW disk gas 
Gaussians and the left-over Gaussians are shown in Fig- 
ure |3t. 

4.3. General Features Observed in the HVG and IVG 
Distributions 

Though a thorough analysis of Figure [7| is beyond the 
scope of this paper, several general characteristics of this 
figure are worth pointing out, especially in consideration 
of similar all-sky maps produced earlier. We spe cifically 
compare our figure to Figure la of |Wakker ( 2004| . Many 
of the same general features may be seen in both maps, 
but some structures are not consistently seen because 
different schemes were used to discriminate MW disk 
from other gas. Wakker employed a symmetric, analyti- 
cal model for the Milky Way disk, whereas we define and 
remove the Milky Way gas empirically. Morever, because 
of the two to four times higher net spatial resolution of 
the LAB data, our Figure [7] map shows some finer struc- 
tures, and our colorscale levels reveal more detail because 
of the > 10 times higher velocity resolution. Most ger- 
mane to the present discussion is that the tendril-like 
structure of the MS and Leading Arm Feature (LAF) 
is more obvious in Figure [7| than in the Wakker maps. 
The same kind of fine-structure is also seen in other HVG 
complexes in our maps (e.g.. Complexes A, C, and M). In 
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Fig. 7. — Sky distribution of IVC and HVC H I gas in Galactic coordinates, where the color represents Vlsr in km s"-*^ and the individual 
Gaussian components are collapsed into single points (compare to Wakker 2004). The Magellanic Clouds (lower left) and Stream (arcing 
across the bottom) as well as Complex C (upper right) and other HVC complexes are clearly visible. 



this regard, our data are more similar to those of P03 and 
B05, with only slightly lower spatial resolution, equiva- 
lent or higher velocity resolution, but covering the entire 
sky. It is because of the higher velocity resolution data 
and the analysis techniques we use to analyze them 
that we are able to build on P03's previous work on tne 
MS. The B05 data, which are higher in spatial resolu- 
tion than the LAB data, were unavailable to us when we 
began this analysis. 

5. RESULTS OF THE GAUSSIAN DECOMPOSITION 

5.1. Magellanic Stream Coordinate System 

Because the Magellanic Stream consistently follows 
such a long trace across the sky, it is useful to have a 
coordinate system for which the equator lies along the 
great circle of the MS. A Magellan ic coordinate sys- 
tem was defined by (Wakker (2001), in terms of the 
great circle along the / = 90 and / = 270° Galac- 
tic meridian. Although this is close to the MS, the 
equator of that coordinate system is not exactly along 
the Stream. We define here a new coordinate system 
which we call the "Magellanic Stream" coordinate sys- 
tem and whose equator more closely bisects the Mag- 
ellanic Stream. The equator of the system was set by 
finding the great circle best-fitting the MS. The pole of 
this great circle is at = (188.5°, -7.5°); the longi- 
tude scale is defined in such a way that the center of the 
LMC b = 280.47°, -32.75°: van der Marel et al. 2002) 
has Lms = 0°. Like the Magellanic coordinate system of 
Wakker, Lms decreases along the MS (towards its tail). 
Many subsequent figures in this paper will use the Mag- 
ellanic Stream coordinate system (Lms, ^ms)- 



5.2. Representations of the Magellanic Stream 

The results of the Gaussian decomposition, with the 
MW disk gas and local zero- velocity gas removed, can 
be seen in Figures |8] and [9) which show integrated inten- 
sity of the Gaussians in three perspectives {Bus vs. Lms, 
Vlsr vs. Lms, Vlsr vs. Bus)- Figure [sji also shows the 
Gaussians of the Magellanic System in Bus vs. Lms with 
color representing <Flsr>. The large velocity gradi- 
ent is evident in this figure. Figure similarly shows 
the sum of Gaussian centers of the Magellanic System in 
^LSR vs. Lms with color representing <Bms>- Some 
readily apparent, prominent features in these figures, 
such as the LMC, SMC, the MS, and the LAF are la- 
beled in Figure [Sj 

5.3. The Two Filaments of the Magellanic Stream and 
their Source 

The two-filament structure of the MS, previously 
pointed out by Cohen (1982) and by Morras (1983), and 
studied recently by P03, is clearly visible and separated 
in our datacube of Gaussian centers. The MS filaments 
can be distinguished between —40° < Lms < 0° iii the 
"^LSR vs. Lms plot (Fig.|9]3). For Lms < -45° it becomes 
difficult to disentangle the two MS filaments in any pro- 
jection of the data and we therefore focus here on the 
^MS > —45° region of the Magellanic System. Figure 
[9] allows us to distinguish both filaments in the head of 
the Stream. Until now the Stream was only known to 
be spatially bifurcated and since the Magellanic Bridge 
gas overlaps the MS filaments near the Clouds it was 
not possible to distinguish or separate the filaments for 
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Fig. 8. — The Magellanic Clouds and Stream H I Gaussians as distributed on the sky. (a) Hue indicates <Vlsr> and intensity indicates 
^Hl (on a logarithmic scale), (b) The three Leading Arm complexes I-III (including the three "clumps" of LA I clearly seen in panel 
(a)), the Large Magellanic Cloud (LMC), Small Magellanic Cloud (SMC), and the Magellanic Stream (MS) are shown in this greyscale 
representation of H I column density, Nm, in units of lO-*^^ atoms cm~^. 



^MS ^ —20°. The velocity bifurcation at the head of the 
Stream ahows us to trace the filaments back further to 
their source than previously possible. One striking char- 
acteristic of the filaments is their oscillatin g pa ttern (Fig. 
and [lO|, which is discussed further in gSTsJ 
Figure 1 10^ is a close-up view of Figure ^ and gives a 
clearer picture of the two filaments (traced with red and 
green lines in the lower- left inset). The filaments cross 



near Lms ^ —28° , but a narrower Bus range (—8.0° < 
^MS < —1.0°) shows the continuity of the "green" fila- 
ment at this point (Fig. [IOJd) and allows us to track the 
filament across this longitude. Near Lms ^ —16° the 
"red" filament crosses the SMC/Bridge gas. Another 
Bus range (—4.5° < Bus < 2.0°) reveals the continuity 
of the "red" filament through this region (Fig.IlOfc). Be- 
yond this point the "red" filament connects to tne LMC 
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Fig. 9. — Integrated intensity (sum of Gaussian areas) of the Magellanic Cloud and Stream H I Gaussians (at their central Vlsr velocity), 
(a) Column density, Nm, in units of 10-*^^ atoms cm"-^ (a blowup of Figure [sj)). A spatial bifurcation of the Magellanic Stream into two 
filaments (first pointed out by P03) can be seen for —40° < i^MS ^ —20°. Other bifurcations are seen farther down the Stream, (b) 
^LSR vs. Lms (hue indicates Bms, and brightness indicates integrated intensity along Bms)- This panel shows the two filaments of the 
Magellanic Stream also to be bifurcated in velocity from —40° < Lms < —20°. One of the filaments is discernible all the way to (LmSj 
^lsr) ~ (—5, +247 km s"-*^) where it connects to the LMC. The second filament can be followed only to (Lms 5 ^lsr) ~ ( — 16.5, +220 
km s"-*^). More velocity bifurcations (or multiple splits) are evident at more negative velocities. The two filaments show strong periodic 
patterns for —40° < Lms ^ after which they follow a fairly linear negative velocity gradient (—95 < Lms ^ —40°). (c) Integrated 

intensity of H I Gaussians in Vlsr vs. ^ms (colorscale in units of K km s"-*^). The bifurcation of the MS is discernible as well as some 
spiraling patterns of the filaments. 
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(hereafter the "red" filament wih be cahed the "LMC" 
filament and the "^reen" filament the "second" filament). 



Putman et al. ( 1998 ) pointed out an "emission filament' 
emanating from the LMC (seen in the sky channel maps), 



but claimed it went into the Bridge; McG ee fc Newton 
|l986) earlier also saw possible indications of an LMC 
filament. However, the P-V representation of Figure [To] 
may be the first conclusive evidence that any part of 
the MS comes from the LMC. P03 claimed that the two 
MS filaments came from the SMC and Bridge, and most 
subsequent tidal simulations (Connors et al. 2006, and 
references therein) have adopted these origins as their 
starting point (e.g., by modeling the SMC as an N-body 
and the LMC as a rigid potential). We argue here that 
the assumptions that are the foundation of these simula- 
tions need to be reconsidered. 

In order to track the LMC filament back to the location 
of its origin in the LMC, we isolated the filament in the 
^LSR vs. Lms plot by eye with the dashed lines shown 
in Figure [TT^. Figure [TTJd shows the distribution of the 
LMC filament on the sky. The filament appears to have a 
spatially periodic pattern. The filament emanates from a 
region of dense H I on the southeastern, or leading, edge 
of the LMC, namely the SEHO (throughout this paper 
we will use the term SE H I overdensity to mean the 
entire region of high-density H I in the southeast of the 
LMC, 05h34"^< a < 05^52"^, -68°28'< S < -71°530. 
The connection to the LMC can be even more clearly 
seen by overlaying the high spatial resolution H I data 
of the LMC from Staveley-Smith et al. (2003; hereafter 
S03) on our map (Fig. 112b). The LMC filament is clearly 
arm "B" seen by S03. 

The SEHO is a natural place for an H I stream to 
originate, due to the high-density of H I there (the high- 
est concentration of H I in the LMC). Furthermore, the 



SEHO is near the end of the LMC bar, is rich in CO (Ya- 
[maguchi et alT||20QT ), Ha emission jKim et aL||l999 ), g: 



ant molecula r clouds ([Yamaguchi et al.|200l| ), and young 
star clusters ( Bica et al.|1999[). T here are also several su- 
pergiant shells plim et al" Tggg) and two CO filaments 
lYamag uchi et al.| |2001|) in this region. Their relation to 
the MS is further explored in ^ 
The origin of the second filament is not as well defined. 



In our maps (Figs.[9]and 10) the second filament can only 
be clearly traced tonigheFlongitude as far as the Magel- 
lanic Bridge (near Lms ^ —15°, V^sr ~ +200 km s~^). 
It is not clear whether the second filament eventually 
connects to the SMC, to the Bridge, or to the LMC. We 
used a by-eye selection in the Vlsr vs. Lms plot (similar 
to the one used for the LMC filament above) to extract 
the second filament gas from the database of Gaussians 
(the dashed lines in Fig.jlTJ^); its distribution on the sky 
is shown in Fig. pT] i. 

Figure [l3^ compares the positions of the two filaments 
on the sky (after the velocity selections in Fig. [TT^ and 
c have been applied). The filaments can be more easily 
distinguished in this representation because the Bridge 
gas has been removed. The two velocity filters (and 
the filaments themselves) overlap in Flsr and Lms near 
^MS ~ — 17° and Lms ~ —28° (as indicated by the dot- 
ted lines). Therefore, in these regions the separation of 
the filaments using the velocity selection alone is not as 
good, but, fortunately, the filaments are spatially well- 
separated in these regions. The continuity of the two fila- 
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Fig. 10. — (a) Close up of the Vlsr vs. Lms distribution of 
Magellanic Clouds and Stream H I Gaussians from Figure lOp (hue 
indicates Bms, and brightness indicates integrated intensity along 
^Ms)- The bifurcation of the two Magellanic Stream filaments is 
clearly discernable as are strong periodic patterns in the radial ve- 
locities of the filaments. The inset illustrates the shape of the two 
filaments (red - LMC filament, green - second filament), (b) A 
narrower range in ^ms to show that the second filament is contin- 
uous at the point (near Lms ~ —28°) where the two filaments cross 
in this projection (on a linear scale from 0-100 K km s"-*^ [white 
to black]), (c) A slightly different range in Bms to indicate that 
the LMC filament is continuous as it crosses the SMC/Bridge gas 
near Lms ~ —16° (on a linear scale from 0-200 K km s~-^ [white 
to black]) 
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Fig. 11. — The two Magellanic Stream filaments isolated by 
velocity. The LMC filament: (a) Vlsr vs Lms distribution 
for the Magellanic Cloud and Stream H I Gaussians showing two 
Magellanic Stream filaments (same as Fig. [ob). The gray dashed 
lines show the velocity limits used to isolate The LMC filament. (6) 
Sky distribution of the column density, A^hi, of the H I Gaussians 
for the LMC filament selected by the velocity limits shown in panel 
(a). The association of the LMC filament with the LMC and the 
spatial periodic patterns are apparent. The second filament: (c) 
Vlsr vs. Lms for the Magellanic Cloud and Stream H I Gaussians 
showing the two filaments (same as Figure |9|)). The gray dashed 
lines show the velocity limits used to isolate the second filament. 
(d) As in panel (6), but for the second filament selected by the 
velocity limits shown in panel (c). The second filament can only 
be distinguished for Vlsr ^ —17° and its source remains unclear. 
The top color bar indicates ^ms for panels (a) and (c), while the 
bottom color bar indicates column density, A^hi (in units of 10-*^^ 
atoms cm~^), for panels (b) and (d). 

ments in these overlap regions can also be seen in Figure 
[Tq] 3 and c. The patterns of the two filaments on the sky 
are similar (also see Fig.|8]3), shifted by in Lms and 
~1° in Bus- Possible p hysic al explanations for these 
patterns are discussed in §5.5| Figure [13]3 shows the H I 
column density on the sky of the Magellanic Gaussians 
with the "guiding" centers of the LMC filament (red), 
the second filament (green) , and the connection of LA I 
to the SEHO (orange, discussed in ^5.4) overplotted in 



order to indicate their positions relative to more familiar 
structures (e.g., SMC, Bridge, etc.). 

5.4. The Source of the Leading Arm Feature (LAF) 

The LAF consists of three complexes of gas north of 
the MCs. Pieces of these s tructures were f i rst se en by 
van Kuilenburg (1972) and [Wannier etlH] (|1972D. Al- 



though there was no direct connection, 'Ma thewson et al.| 
(Il974) suggested a possible association of these features 
to the MCs and th is hypothesis w as furt her e xplored by 
Mathews on et al.| ([Ty9), Morras (1982"), and 'Ba.i a.ia"et 



al. ( 1989 ) . With reprocessed HIPASS data Putman et al. 
(jl998) demonstrated that the LAF gas is indeed associ- 
ated with the Magellanic Clouds and is an extension of 
the Magellanic Stream. The association of the LAF with 
the Magellanic Clouds and the MS lent support to the 
tidal origin of the MS over the ram-pressure model, be- 
cause a leading feature of the MS would not be expected 
if ram-pressure were the dominant force. 
The three primary complexes of the LAF can be seen 
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Fig. 12. — Close-up of the Fig. Ill| 3 map of the integrated in- 
tensity (sum of Gaussian areas) oi the LMC and LMC filam ent 
H I Gaussians on the sky after a velocity filter (see Fig. |ll^ ) is 
applied. These maps show that the LMC filament is emanating 
from the SEHO in the LMC when viewed with either (a) the LAB 
data only, or {h) the high-resolution H I data (255.0 < Vlsr < 
288.0 km s"-*^) from S03 are substituted in the region outlined by 
the dashed lined box. The filament can be associated with S03 
arm B. A square root transfer function is used in these greyscale 
images. 

in Figure Is] at positive Lms: LA L (3° < Lms < 29°, 
-34° < 5ms < -6°); LA IL (36° < Lms < 61°, -17° < 
^MS < -10°); and LA IIL (35° < Lms < 62°, -2° < 
^MS < 11°) (nomenclature by B05). Our analysis here 
focuses on LA I, the closest LAF complex to the Magel- 
lanic Clouds. LA I consists of three nearly rectangu- 
lar ~2.5° X 8° concentrations of gas that each lie al- 
most parallel to lines of constant Lms, and, combined, 
form a linear structure ma king an ^40° an gle with the 
^MS=0° line (Fig.|8|. Put man et ar| ( |1998| ) showed that 
the two concentrations nearest the LMC are nearly con- 
tinuous (see their Fig. 3) and it is therefore likely that 
the entire LA I feature is a physically connected struc- 
ture. The first concentration of LA I is close to the 
south-eastern edg e of the LMC , both in position and 
in radial velocity. Putman et al. (1998) claim, however, 
that the Leading Arm material comes mainly from the 
SMC, based on a filamentary feature that is nearly paral- 
lel with the 5ms=0° line and that begins near the SMC 
and stretches to the first concentration of LA I (see their 
Fig. 1). However, S03 noted several H I features of the 
LMC (arms B, E, S, and W) and remarked that arm E 
pointed to the Leading Arm clouds, which lay beyond 
the coverage of their survey. S03 go on to say that deep, 
reprocessed HIPASS data (P03) shows a continuous con- 
nection between arm E and the Leading Arm. Never- 
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Fig. 13 . — (a) The sky distribution of the two Magellanic Stream filaments as we have extracted them with velocity filters shown in 
Figure [TTk,c (using a logarithmic transfer function for the intensity). The LMC filament is shown in red and the second filament in green. 
Most ot tne LMC is included in the velocity filter of the LMC filament. The two velocity filters (and the filaments) overlap in Vlsr and 
Lms near Lms ~ — 17° and Lms ~ —28° (as indicated by the dotted lines). Therefore, in these regions the separation of the filaments 
using the velocity selection alone is not as good. However, in these regions the filaments are spatially well-separated (the continuity of the 
two filaments in these overlap regions can also be seen in Figs. |10| 3 and c). (6) Column density, A^hi, of the H I Gaussians (on a logarithmic 
scale from 1-200x10-*^^ ato ms c m~^ [white to black]). The "guiding" centers of the LMC filament (red), the second filament (green), and 
LA I (orange, discussed in ^5.4| are overplotted in order to indicate their positions relative to more familiar structures (e.g., SMC, Bridge, 
etc.). 

theless, S03 conclude that the Leading Arm gas mainly 
arises from the SMC, and only some LMC gas "leaks" 
into the Leading Arm. B05 showed that the first con- 
centration of LA I is directly connected in position and 
velocity to the H I clouds close to the LMC, but claim 
that it is associated with the Bridge. We believe that 
there is a much firmer association of the Leading Arm 
with the LMC. 

The radial velocities of the LA I complexes are quite 



similar to those of the LMC (Fig.[4]3), but typically some 
~100 km s~^ higher than those oithe SMC. This can be 
seen especially well in Figure [MJd, where the first con- 
centration of LA I (-17° < Bus ^ -10°) is connected 
in position and velocity to an extension of the LMC 
(-10° < Bus < -5°), namely SOS's arm E. To better 
probe the association of LA I to the LMC we show the 
distribution of gas in the velocity range 250 < Vlsr < 
320 km s~^ in Figure 15 i. In Figure [l5p we also overlay 
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Fig. 14. — Integrated intensity distributions of the Lead- 
ing Arm Feature (LA I) and some LMC H I Gaussians: (a) 
^LSR vs. Lms distribution (hue indicates <Bms>^ brightness indi- 
cates integrated intensity along ^ms): (^) ^LSR vs. ^ms (hue 
indicates <Lms>, and brightness indicates integrated intensity 
along Lms)- The thr ee c oncentrations and the LMC arm E are 
labeled (also see Fig. 15fc. These two fi gures show that the first 
concentration of LA I (—17° < ^ms ^ —10°) connects in position 
and radial velocity to an extension of the eastern part of the LMC 
(-10° < Bms < -5°), which is arm E of Staveley-Smith et al.| 
( |2003i .^ ^ ' 

the 803 high resolution ATCA H I data over our own (us- 
ing the same velocity cut) to confirm the basic geometry 
seen in the LAB data, namely SOS's arm E extending out 
of the LMC and towards the first two concentrations of 
LA I (which are not covered by the 803 survey). Arm 
E has the same elongated shape (parallel to Bus) and 
continues the spatial progression of the three LA concen- 
trations (more negative ^ms towards higher Lms)- Arm 
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Fig. 15. — Maps of integrated intensity (sum of Gaussian areas) 
of the Leading Arm Feature (LA I) and the LMC H I Gaussians on 
the sky for gas with 250 < Vlsr < 320 km s"-*^ (on a logarithmic 
scale): (a) LAB data only; (b) LA B data augrn ented with the high- 
resolution H I data from Stavele y-Smith et al.| ([2003) in the dashed 
lined box. The three concentrations ot LA 1 and the LMC arm E 
are marked. These maps show that the first concentration of LA I is 
close to the southern end of LMC arm E. Arm E (which originates 
in the SEHO) is similar in shape to the three concentrations of LA 
I and continues the spatial progression. This indicates that LA I 
originates in the SEHO of the LMC. 

E also continues the velocity trend with Bus as seen in 
Figure [mJd. There is a gap of a few degrees between 
the end of arm E and the beginning of the first concen- 
tration of LA I (although there is a small clump of gas 
between them at [Lms, ^ms]^[4°, — 11°]). However, our 
data also show gaps between the three concentrations of 
LA I that the deeper HIPASS data show are contiguous. 
Therefore, it is likely that arm E and the first concentra- 
tion of LA I are also connected. For all of these reasons, 
we strongly suspect that LA I is physically connected to 
arm E (which starts in the SEHO) and has its origins in 
the LMC. Therefore, we conclude that both the trailing 
LMC filament of the MS and LA I have their origin in 
the SEHO of the LMC. We discuss the implications of 
these findings further in ^ 
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As previously mentioned, IPutman et al. (1998) argue 
that the LAF comes from the SMC, based on a horizon- 
tal feature south of the LMC which seems to connect the 
SMC to the Leading Arm (see their Figs. 1 and 2). This 
feature is extended nearly parallel to the 5ms =0° line 
and the first of these features is also visible in our maps 
(see Fig.[9fi; -10° < Lms < 1° and -8° < Bus < -T) 
and does connect to the end of arm E in position as well 
as in velocity. It is this apparent connection, as well as 
the presence of the other "horizontal" cloud between the 
first feature and the SMC, that is the basis for Putman 
et al.'s claim that the LAF originates in the SMC. We 
cannot definitively refute this hypothesis but believe that 
the geometry of the gas in LA I and arm E is more sugges- 
tive of the LMC as the originator of the LAF: whereas 
arm E is already oriented in the same direction as the 
LA I concentrations and continues their staggered verti- 
cal striping spatial pattern (Figs, [s] and 15), the Putman 
et al. horizontal features are oriented ortnogonally. 

If this gas is moving towards the LMC, then it might 
be possible that the SMC is contributing some gas to 
the Leading Arm or creating its own leading arm. On 
the other hand, if the horizontal feature is moving from 
the LMC to the direction of the SMC, then it may be 
originating from the end of arm E and is part of the 
trailing MS. The structures and patterns of the horizon- 
tal features look much like the filaments of the MS, and, 
when extrapolated, the horizontal features seem to con- 
nect to the second MS filament (Fig. |9^). If this is the 
case, then the entire Magellanic Stream and the Leading 
Arm might originate from LMC gas; more work is needed 
to clarify this part of the Magellanic Bridge. Based on 
the evidence that we have presented in this section (sum- 
marized in Figs. p!4| and p!5|), however, we conclude that 



most (if not all) of the Leading Arm gas originates in the 
SEHO of the LMC. 

5.5. The Cause of the Periodic Pattern 



As pointed out in §5.3[ the two filaments of the MS 
have pronounced sinusoidal patterns in velocity and in 
position, especially for Lms > —40° (Fig. |9|. P03 noted 
that the two filaments give the impression of a double 
helix and postulated that it might be due to the "pseu- 
dobinary motion of the LMC and SMC" . This tumbling 
motion thus provides a first hypothesis as to the cause 
of these patterns. However, there are several problems 
with this hypothesis. The path of two bodies orbit- 
ing each other would create a double-helix pattern in 
our position-position-velocity (ppv) datacube and in the 
projections of the datacube. That is not seen for the 
two filaments at the head of the stream (Lms > —40°), 
although they might cross farther downstream. Figure 
16 shows the paths of the central concentrations of the 
two filaments, in the three projections (5ms vs. Lms, 
"^LSR vs. Lms7 and li,sR vs. 5ms) ^ and reveals that the 
two filaments do not cross or spiral around each other 
(at least in this Lms range), but rather create two in- 
dependent spirals parallel to each other. Therefore, the 
spiraling motion that we see cannot be explained by the 
tumbling of the LMC and SMC about each other (if 
they are indeed bound to one another). It is also not 
clear where the second filament originates, whereas its 
origin must lie in the SMC if the tumbling hypothesis 
is to work. The two filaments do eventually cross each 




Fig. 



16. — This figure is analogous to Figure [9] The path of the 

irom three perspec- 



central concentrations of the two MS filaments 
tives: (a) Bms vs. Lms, (^) ^lsr vs. Lms, and (c) Vlsr vs. Lms- 
The LMC filament (solid) and the second filament (grey, dash- 
dotted) do not wrap around each other, but move parallel to one 
another, in this Lms range. 

other at Lms ^ —47° (Fig. [9|i) and other places, beyond 
the large velocity oscillations seen at the head of the MS 
(-40° <Lms < -5°). The spatial crossing of MS strands 
we see at Lms ^ —47° has become apparent in our maps 
because the zero-velocity gas has been removed; this MS 
crossing is also apparent in Figure 2 of B05, but was not 
noted by P03. It is possible that these "later", more 
widely separated from the MCs, crossings could be due 
to the binary motion of the Magellanic Clouds (if the 
second filament actually originates in the SMC), but the 
spirals at the head of the Stream cannot. There must be 
another explanation for these spiraling patterns close to 
the MCs. 

We postulate instead that the spiraling motion is an 
imprint of the rotation curve of the LMC (and possi- 
bly the SMC). This scenario fits the data better than 
the previous hypothesis, especially for the LMC filament 
where the sinusoidal variation is pronounced and the fil- 
ament can be reliably traced back to its oriffln in the 
LMC, which has a well-defined disk (e.g., Weinberg & 
Nikolaev 2001; van der Marel 2001) and rotation curve 



(e.g., |Kim eFal ||1998| IvaTHer'T^ al.|20Q2| . The 
amplitude of the spatial variations is ^2 (Fig. fTTl)) 
which is close to the radius of the SEHO from the LMC 
center. A sinusoid+line fit to the LMC filament in 
lisR vs. Lms gives 

Flsr= (26.4kms-i) sin(f^LMS-5.9°) + (2) 
(4.40 km s-^ deg-^) Lms + 291.75 km 

and can be seen in Figure [iTj Assuming that the linear 
portion is due to the orbital motion of the LMC about 
the MW, the sinusoidal amplitude of 26.4 km s~^ is not 
that different from the projected LMC rotation velocity 
at the radius of the SEHO (~2.5° from the center of the 
LMC), i.e. Kot sin i ^ 39.4 km s"^ sin 34.7° = 22.4 
km s-^ (vdM02; |Kallivayalil et al. 112006811). It therefore 



seems that the data are reasonably consistent with the 
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Fig. 17. — Vt. s-r vs. Lms of the Magellanic Stream Gaussians 
(same as Fig. |10| with the sine+linear fit (Equation 2) to the LMC 
filament. 

hypothesis that the spiral pattern in the LMC filament 
was created by the rotation of the SEHO (the filament's 
apparent birthplace) about the center of the LMC as the 
entire system orbited the MW. Our hypothesis is that 
some force — tidal, ram pressure or otherwise — must 
have pulled/pushed the H I gas out of the SEHO (see 
Q and the trailing gas bears the imprint of the SEHO's 
particular velocity and position. 

Assuming this hypothesis is true, we can estimate the 
drift rate of the H I gas in the LMC filament by us- 
ing the angular period of the sinusoidal pattern. The 
orbital period of the SEHO around the LMC is T = 
(27r X 2.5°x0.875 kpc deg-i)/39.4 km s'^ ^ 340 Myr. 
In this time the SEHO's radial velocity undergoes one 
full cycle, which we can identify with one full cycle of 
the LMC filament's velocity pattern. Even though we 
have a spatial period from the sinusoidal fit, it is ap- 
parent from Figure [T7| that the sinusoid+line deviates 
slightly from the periodic pattern. We can get a more 
accurate estimate of the spatial period by looking at the 
maxima at (^^ms^^lsr) = (—10.3°, 271.4 km s~^) and 
at (Lms^lsr) = (-29.8°, 194.6 km s-^), which gives 
19.5° for the spatial period. The LMC filament gas is 
probably still close to the LMC distance from the Sun 
of ~50 kpc which we can use to calculate the distance 
between the maxima as I) = 19.5° x 0.875 kpc deg~^ = 
17.06 kpc. If the time difference between when these two 
points of gas were in the SEHO is just 340 Myr, the or- 
bital period of the SEHO, and if the drift rate of the LMC 
filament gas away from the LMC is also roughly con- 
stant, then we can estimate it as Fdrift = 17.06 kpc/340 
Myr = 49.1 km s~^. This is a lower limit since it does 
not include any radial component of motion. The drift 
velocity is ~l/9 the tangenti al velocity of the LMC, 
'^LMC,tan = 367 km s~^ ([Kallivayalil et al. 2006a). 

Lastly, we can use the drift rate of the LMC filament 
gas to roughly estimate the age of the Magellanic Stream. 



If we assume that the drift rate has been approximately 
constant for the entire MS as well as an average distance 
of 50 kpc to the Stream, then the age of the Stream is 
(100° X 0.875 kpc deg-i) / 49.1 km s'^ = 1.74 Gyr. A 
more direct way of calculating this would be to count the 
number of oscillations along the MS and multiply by 340 
Myr (the orbital period at the SEHO radius). However, 
the amplitude of the modulations decreases quite dra- 
matically after about ~40° from the LMC and the two 
filaments are not as well separated in velocity so it is not 
yet possible to do this in practice. If the distance to the 
tip of the MS is larger than 50 kpc (as predict ed by the 
models of [Connors, Kawata, &: Gibson] [20061 see their 
Fig. 6), then the MS would be older, and conversely, the 
MS would be younger if the MS tip is closer tha n 50 kpc 



(as suggest ed by the ram pressure models of Mastropietro 
et al.|20Q5[ see their Fig. 7). A changing drift rate would 



also alfect the MS age. Even so, an MS age of ~1.7 Gyr 
is fairly close to the 1.5 G yr age of the Stream found by 
Murai & Fujimoto' ("1980) at which time their tidal sim- 



ulations show an SMC-LMC-MW close encounter. How- 
the orbit calculations using the new HST proper 



ever, tne orbit calculations usmg 
motions f Kallivayalil et al.|2006a b 
szewski 20071) and an JNFW MW pc 



Piatek, Pryor, & 01 



esla et al. 



7^ and an JNl^^W MW potential by^ 

(2007) give past MC encounters at roughly (3.2, 3 and 6 
Gyr ago (G. Belsa, private communication). These num- 
bers are highly dependent on the mass ratio of the MCs 
(~1:10) and the SMC's proper motion both of which are 
quite uncertain. A close encounter at ~2-2.5 Gyr would 
be consistent with the bursts in the star formation rate 
seen in the star formation histories of both MCs ( S meeker-] 



Hane et al j2002[jHarris fc Zaritsky|2004| ). Therefore, an 



age of ~ 1.7 Gyr for the MS seems reasonable considering 
the uncertainties. 

Even though the rotation hypothesis seems self- 
consistent there is one possible problem - the coherence 
of the sinusoidal pattern in the filaments. If the rota- 
tion hypothesis is correct then there should also be an 
oscillating pattern in the transverse velocity of the fil- 
aments. One would expect that the variation in veloc- 
ities would pull apart or stretch the filament. At the 
very least the different velocities should place the gas 
in different MW orbits because of their differing ener- 
gies, which would also destroy the coherence of the fila- 
ment. Since this is not seen, and the structure is coherent 
over at least ~40°, there mus t be something else hold- 
ing the structure together. |Konz, Briins, fc Birk (|2002| 
show how the interaction of a cold cloud with its sur- 
rounding ambient medium can create a magnetic barrier 
which helps keep the cloud intact. In addition to explain- 
ing the survival times of the MS filaments, the interac- 
tion of the clouds with their surrounding medium might 
help explain the coherence of the Stream over such large 
distances. More detailed N-body/hydrodynamic (and 
magneto-hydrodynamic) simulations are needed to inves- 
tigate how the MS filaments, with their oscillating veloc- 
ity patterns, can remain coherent over such distances. 

For Lms ^ —40° the sinsusoidal pattern diminishes 
substantially, and the amplitude of any leftover varia- 
tions is much smaller (Fig. [9| and the behavior of the 
MS is more nearly linear {Vlsr ~ 7.60 Lms +370 km s~^; 
for Lms < —40°). Something dramatic must happen at 
^MS ^ —40° to cause this change. One possible expla- 
nation is that this is where drag from the MW halo gas 
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Fig. 18. — Integrated intensity of the LMC H I datacube (from 
803) from three perspectives: (a) Column density, A^hi? in units 
of 10-*^^ atoms cm~^, (b) Vlsr vs. a (integrated in S; greyscale in 
units of K), and (c) Vlsr vs. 6 (integrated in a; greyscale in units 
of K). From (c) it is clear that there is gas coming off the LMC 
on the high-velocity end and moving to the south. The CO clouds 
from Mizuno et al. (2001) (white dots in all panels) are almost 
always associated with the LMC disk component (the regions with 
higher integrated intensity in panel (c)). 



becomes important, dampens the sinusoidal pattern, and 
causes the MS to fohow a quite hnear curve. Another 
possibihty is that the LMC is actuahy much larger than 
previously thought (exceeding even the newly suggested 
limit in Muhoz et al. 2006) and the MS gas is escaping 
the tidal radius of the LMC at Lms ~ -40°. 

If the LMC tidal radius is 40 kpc then from equation 
(45) in vdM02 (and a MW mass interior to the LMC of 
4.9 X 10^^ Mq) this would imply a minimum LMC mass 
of L68 X 10^1 Mo (or 34% of Mmw(< 50 kpc)). The 
total intrinsic visible luminosity of the LMC is Ly = 
3.0 X 10^ Lq (vdM02) which implies a mass-to-light ratio 
of M/Ly = 56. This would mean that the LMC is very 
dark matter dominated. 

If we assume that ram pressure is the dominant force 
causing the MS to drift back behind the LMC then, 
with some other assumptions, we can estimate the den- 
sity of the hot MW halo gas. The ram pressure on 
the MS is P = Pmw^ms- -"-^ ^ms is the approxi- 
mate diameter of the MS, then the acceleration that 
the MS experiences is aus ^ Pmw^ms^ms/(Pms^ms) ^ 
Pmw'^ms/(/^ms^ms)- Solving for the ratio of densities 
gives Pmw/pms ^ <^ms^ms/^ms- assume that the 

MS was undergoing a constant deceleration due to ram 
pressure in the recent past then we can estimate this 
deceleration from the sinusoidal pattern in the LMC fil- 
ament, aus ~ 2Ax/At2 ^2x 17.06 kpc/(340 Myr)^ ^ 
295.2 kpc/Gyr^. Approximating the diameter of the 
LMC filament at the head of the Stream as Lms ^2 
kpc and ^ms ^ '^lmc = 378 km/s (|Kallivayalil et 



al.| |2QQ6a 
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LMC filament is A^hi ~ 1 x 10^^ atoms/cm . Assum- 
ing a distance of 50 kpc and a width of ~2 kpc gives 
a number density of nus ^ 0.016 atoms/cm^. Finally, 
we derive the number density of the hot MW halo as 
nMw ^ 6.3 X 10~^ atoms/cm^. This rough estimate of 
the density of the hot MW h alo gas is consistent with 



most previous estim ates (e.g. Stanimirovic et al. 2002 
Sembach et aL]|2003D . 



6. THE CAUSE OF THE OUTFLOW FROM THE SE H I 
OVERDENSITY 

6.1. LMC High Velocity Gas Ejection 

What is causing the Leading Arm and LMC filament 
gas to flow out of the SEHO? It is difficult to assess 
this question with the LAB data due to its relatively low 
spatial resolution (0.5°); however, the higher spatial and 
velocity resolution of the S03 H I data of the LMC from 
the Parkes telescope is well suited to a closer look at 
the SEHO. In Figure 18 we show the integrated intensity 
of the LMC H I from three different perspectives: (a) 
5 vs. a, (b) Vi^sR vs. a, and (c) Vlsr vs. S. The CO 
clouds identified with NANTEN by [Mizuno et al.| ( |2001 ) 
are plotted as white dots in these figures. 

The high density of the SEHO is readily apparent in 
the column density plot (Fig. IlSk), as are great voids in 
the gas, which are evidence oTsupergiant shells.^ The 
Leading Arm can be seen stretching southward of the 

^ These supergiant shells are even more pronounced in the much 
higher spatial resolution ATCA data used by Kim et al. (1999) 
in their study of the supergiant shells (as well as in the combined 
ATCA+Parkes H I data by Kim et al. 2003). 



18 



Origin of the Magellanic Stream 



200 400 600 80 250 500 750 1000 1250 15 00 

350^ ■ ■ ■ ■ ■ 
(b) 04:27:52 < RA < 05:32:52 ^ 

Center/West ^ '' 




-76 -74 -72 -70 -68 -66 
Decllnotlon 



-76 -74 -72 -70 -68 -66 
Decllnotlon 



Fig. 19. — Integrated intensity of the eastern and central/ western 
regions of the LMC from the S03 datacube. (a) The eastern re- 
gion of the LMC showing the high and low velocity components 
(for —72°<6< —69°). (b) The central and western regions of the 
LMC showing a velocity distribution with less outflow. The high- 
velocity gas is coming mainly from the eastern part of the LMC. 
The grescales are in units of K. It is clear that the lower- velocity 
component in the south-east is the "disk" component by its high 
integrated intensity (compared to the high-velocity gas) and its 
similarity to the velocity distribution to the west. The velocity 
criterion used to separate the hi gh-v elocity gas from the LMC disk 
gas at each declination (see Fig. |20| is shown by the dashed line. 

SEHO and the LMC filament diagonally to the south- 
west (SOS's arms E and B, respectively). A position- 
velocity diagram of the LMC (Fig. JTsb; integrated along 
S (which is nearly parallel to the LMC H I kinematical 
line-of-nodes,^ Kim et al. 1998) shows that the LMC is 
not symmetrical in H I but is lopsided and contorted. 
Furthermore, there is a large build-up of neutral gas on 
its leading/eastern ed ge (where the S EHO is located). It 
has been suggested by de Boer et al.| ([l998 ) that this is a 
result of the LMC interacting with the diff use M W halo 
gas (see further discussion of this point in ^7.3). 

The velocity profile of the LMC in the third view (Fig. 
[iS]^), a position- velocity diagram integrated along a (i.e., 
nearly perpendicular to the Kim et al. line-of-nodes) ap- 
pears, overall, to be much more symmetrical because it 
reflects very nearly the LMC rotation curve. However, 
one can see in Fig. 18 3 some gas at higher velocities than 
the LMC disk that has no corresponding, lower velocity 
counterpart. This gas appears to come off of the LMC 
disk (at the high- velocity side) and to stay at higher ve- 
locities than the disk. Both the high- velocity gas and the 
LMC disk gas follow a negative velocity trend with S (i.e. 
lower velocities toward the south), but the high- velocity 
gas remains at higher velocities than the LMC disk at all 
positions. 

In Fig. [19] we compare the eastern versus the central- 
western portions of the LMC (where we have integrated 
along a). The first panel reveals the "high- velocity" gas 
clearly, with most of it coming from the south-eastern 
portion of the LMC, especially from a > 05^30"^ and 5 
< —70°. We isolate the high- velocity gas by making a 
rough velocity selection (the dashed line in Fig. 19) in 
the VhSR vs. 5 plane, and show its sky distribution in 
Figure [2Ql3. This high- velocity gas is indeed the Leading 
Arm and LMC filament (the same as S03's arm E and 
B), and it is clearly coming from the south-eastern region 

^ Here we adopt the deflnition (used by Luks &; Rohlfs 1992 
and Kim et al. 1998) of the line-of-nodes as the axis of maximum 
velocity gradient. This differs from the deflnition of the line-of- 
nodes by vdM02 as the intersection of the galaxy plane with the 
plane of the sky. 



of the LMC, or the SEHO. 

Another obvious feature in Figure [19^ is a strong out- 
flow of gas starting at —70° and moving to higher 
velocities and southward; this feature might be coming 
from a supergiant shell (there are several in that regi on). 
The h igh spatial resolution H I data from Kim et al.| 
( |2003D show this region to have many thin outflows from 
the SEHO (see their Fig. 5, especially 05^37"" < a < 
05^48"^) that are not resolved in the S03 data. As Kim 
et al. point out, these are probably outflows from several 
smaller giant shells (they identified 16 giant shells in this 
region). In their figures it appears that there are many 
small fountains of gas coming out of the SEHO which 
are contributing to the Leading Arm and trailing LMC 
filament. These outflows are explored in more depth in 
t he next subsectio n. 

Luks fc Rohlfe] ^1992[ ) found that the H I gas in the 
SEHO is disturbed and that there are two velocity com- 
ponents, which they called the high- velocity disk ("D") 
and low- velocity ("L") components. Figures 18 -20lshow, 



however, that the high- velocity component in the south- 
eastern LMC is not the LMC disk, but an outflow that is 
the beginning of the Leading Arm and the trailing LMC 
filament. The low- velocity component has a much higher 
column density than the high-velocity component (Figs. 
20 i and c), and is contiguous with the rest of the LMC 
^k in position- veloci ty slices (Fig. [2Tt. Also, the CO 
clouds in this region dMizuno et al. pQOT ]) are all asso- 



ciated with the low- velocity compononent and none are 
associated with the high-velocity component (Figs. 18 



and 
Rho 



21|. Therefore, we conclude (contrary to Luks 6? 
Ts) that the low-velocity gas is the disk component 
(albeit a disturbed disk) in the south-eastern part of the 
LMC. 

6.2. Source of the High-Velocity LMC Gas 

What might the source be of this high velocity out- 
flow? It has long been known that there are high- velocity 
gaseous flows emanating from star formation regions in 
the LMC, particularly in the form of supergiant shells 
(e.g., Meaburn & Blades 1980; Meaburn 1984; Meaburn 
et al. 1987); the most recent com prehensive s t udy o f 
these structures has been done by Kim et al. (1999). 
Several of these supergiant shells are also visible in the 
Parkes H I datacube and are, we believe, relevant to the 



origin of the MS. The a =05^43"'04" slice in Figure |2lji 
shows a bubble feature with a void of gas in its cen- 
ter (-70.5°< 5 < -69.0°, 240 km s-^ < Flsr < 280 
km s ^). This void of gas can also be seen in channel 
maps of the LMC ( Fig. [22) -r; e.g., the U-shaped feature 
most obvious in Fig.[22]i)and near the position of the su- 
pergiant shells SGS 19 and 20 (Kim et al.|1999|). In fact 



[.I1999| . 

the bubt 



the void is probably SGS 19 and 20 and the bubble their 
combined envelope. There is a blob of gas that is "at- 
tached" near the edge of this bubble (Figs. 21 i,e and Fig. 



22b-g). The blob and this linking gas has a "tadpole" - 
Eke structure with the head farther south than the tail; 
together these features have an upside-down L shape and 
end near the position of the southern edge of SGS 20 
(Figs. 21 ",g). This seems to be evidence of an outflow 
from SGS 19 and 20. Figure 23] shows three perspectives 
of this outflow. SGS 1 also has an attached blob of gas 
that exhibits the same tadpole features (see Fig. 



24 



^6.4). We conclude that both of these are outflows 



2 and 
3eing 



Nidever et al. 



19 



10 



100 



10 



100 



c 
o 

-I—' 

o 

Q 



-62 
-64 
-66 
-68 
-70 
-72 
-74 





































\ \ \ 


rr 





7.0 6,5 6,0 5,5 5,0 4,5 4,0 
Right Ascension 



c 
o 

-I— ' 

o 

c 

o 

Q 



-62 
-64 
-66 
-68 
-70 
-72 
-74 









































\ \ \ 







7,0 6,5 6,0 5,5 5,0 4,5 4,0 
Right Ascension 



c 
o 

-I—' 

o 

73 
Q 



10 



100 



■62 
■64 
■66 
■68 
■70 
■72 
■74 









gos ~ ' 7^ 






















Ar^^^E — ■ 













7.0 6.5 



6.0 5.5 5.0 4.5 
Right Ascension 



4.0 



c 
o 

-I—' 

o 
c 

"o 

CD 
Q 




7,0 6,5 



6,0 5,5 5,0 4,5 
Right Ascension 



4,0 



Fig. 20. — The sky distribution (column density, Nm, in units of lO-*^^ atoms cm~^) of the different velocity components of the LMC: 
(g) All Magellanic-associated gas, (b) LMC disk gas only, and (c) high- velocity Magellanic gas only, using the velocity selection from Fig. 
|19| Arm E and B are part of the high- velocity component, (d) The distribution of both the LMC disk (red) and high- velocity (green) 
components of the LMC. Near the LMC disk the high- velocity gas is predominantly on the eastern/leading edge of the SEHO. 



expelled by the supernovae (SNe) and high stellar winds 
coming from the supergiant shells. A possible explana- 
tion for the tadpole-like structure is that once the outflow 
reaches a certain height above its origin site other forces, 
such as ram pressure and tidal forces, become dominant 
and push/pull the outflow southward (in the same direc- 
tion that arms B and E are being pushed/pulled). 

We propose the dynamical kick of SNe-driven SGSs as 
the mechanism th at gave rise to th e majority of the out- 
flow in the SEHO. |Meaburn| ( |1984| previously found the 
30 Doradus nebula to be a primary source of outflowing 



gas, but the Parkes data cube reveals this phenomenon 
to be more widespread, encompassing the entire SEHO. 

6.3. Tracing Arms B and E 

The H I arm B from the LMC can be traced back to 
its origin in the region around SGS 18 and 20, as can 
be seen in Figures [22]P-r and Figure [25j There are two 
main strands of arm B, a northern and southern strand 
(Figs. 20 i, c), each very thin (about ~0.25° wide) but 
with distinct velocities. The northern strand appears to 
emanate around (a, 5, Flsr) ~ (05h38"^32% -70°39'16'^ 
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Fig. 21. — Posit ion- velocity cuts in a of the eastern part of the S03 LMC H I datacube (going from east to west). Two a bins in the 
datacube were combined to create each panel. The central ol is shown in the upper-left hand corner of each panel and the greyscale is in 
units of K. The dashed line shows the average rotation curve from the western part of the LMC, which helps separate the high-velocity 
gas from the disk component. White crosses indicate the central position and velocity of the three sugergiant shells from Kim et al.| (|1999| 
that are in this region (SGS 18, 19 and 20; the shells extend to the neighboring panels as well). This figure shows the large outnow of 
gas from this region to high velocity. A large bubble of gas, most clearly seen in panels (c)-(e), (centered at 5=— 69°28^11^^, Vlsr =265 
km s"-*^) is near the position of S GS 19 and an outfiow blob of gas, most apparent in (d)-(g), seems to be coming off the bubble. The CO 
clouds from |Mizuno et aT] ( |200l| (white filled dots) are clearly not associated with the high- velocity outfiow gas but rather with the disk 
component. I he central position and velocity of the 30 Doradus star cluster R136 is indicated by the white box in (g). R136 is a little bit 
west of where most of the high- velocity gas is being blown out. 



273 km s-^) = (5.6422°, -70.6544°, 273 km s-^) (Fig. 



22 
T8 



and Fig. 25 



3, Cj 



near the intersection of SGS 20 and 
The southern strand of arm B is at lower veloc- 
ity than the northern strand (by ~30 km s~^) and ap- 
pears to emanate around (a, ^, Vlsr) ^ (05^38"^09®, 
-71°07W, 200 kms-i) = (5.6358°, -71.1189°, 200 



km s-^) (Fig. [25^) in SGS 18. In Figure [25^ there ap- 
pears to be a second parallel filament at the same ve- 
locity as the southern strand farther to the south. Since 
we don't have distance information there are projection 
effects which make it difficult to pinpoint exact origins 
for the arm B filaments. Therefore, the locations given 
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Fig. 22. — Channel maps of the eastern part of the S03 LMC H I datacube (from high to low velocity). Six velocity channels (^5 
km s"-*^) were combined to form each panel. The central Vlsr velocity (km s"-*^) is shown in the lo wer- left hand cor ner of the panel. Panel 
(a) shows the column density (for 208-291 km s"-*^) with the supergiant shells and giant shells from Kim et al.l(ll999) overplotted in black. 
The scale at the top indicates the column density, A^hi, for panels {b)-{r) in units of 10-*^^ atoms cm^^ White-black dashed lines indicate 



the th ree sugerg i ant sh ells SGS 18, 19 and 20 (at their central velocity, but they extend to many of the neighboring panels in velocity) 
| ( |1999| >. The outlines of these three supergiant shells are also shown as black dashed lines in panels {b)-{r) as reference 
positions. I'h js B gure shows the large outflow of gas from this region to high velocity including arms E and B. A large outflow blob of gas 



from Kim et 



(also see Fig. |21| ), most clearly seen in panels (d)-(g), is near the position where SGS 18 and 20 intersect, but (Mizuno et al. 2001; white 
filled dots in all panels) are located in regions of dense H I and from which the outfiows of gas originate (but at lower velocities than the 
outfiows). The 30 Doradus star cluster R136 is indicated by the white box in h (at its central velocity). 



above should be considered rough estimates. However, it 
is clear that arm B originates in the SEHO region. Arm 
E is wider, more diffuse, and on average at higher veloc- 
ity than arm B. Arm E appears to emanate from many 
outfiows across the SEHO region. There is a trend that, 
moving eastward, the outflow sites move to the north 



and to higher velocity, which is also a general trend in 
the SEHO. 

It is unclear why the two arms are moving in the di- 
rections they are (arm E to the south, and arm B to 
the southwest) even though they originate in the same 
place. The two outflow blobs mentioned above show that 
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the gas must first reach a certain positive velocity offset 
(-30 km s"^ for the SGS 1 outflow, and -20 km s"^ for 
the SGS 20 outflow) from the systematic SGS or LMC 
disk velocity before it starts to move appreciably from 
its place of origin. This might tell us something about 
the forces operating on the gas. The motions of the two 
arms /filaments might be related to which side of the disk 
they are on. Since the LMC disk is inclined relative to its 
direction of motion through the MW halo any gas blown 
out in "front" of the disk (towards the direction of LMC 
bulk motion) will experience more ram pressure than gas 
blown out "behind" the disk (away from the direction of 
motion). All of the outflow gas is at higher positive ve- 
locities and it would seem likely that the gas is on the 
opposite side of the LMC from our perspective and in 
"front" of the LMC as it moves through the MW halo 
(this is discussed further in ^6.7|) . More investigation is 
needed to explain the motions of the two filaments. 

6.4. Mass Accounting 

We can estimate the mass of H I gas in the outflow blob 
(Fig. 23) as well as the mass evacuated from the SGSs 
and see how they compare. The mass of the blob of 
gas (within 05^39"^32^ < a < 05^47"^37% -71°26'02''< 
5 < -70°46'53'^ 267.4 km s-^ < Flsr < 280.6 km s-^) 
is -1.5 X 10^ Mo (using an LMC distance of 50 kpc). We 
can estimate the mass of the H I gas that has been evac- 
uated from SGS 19 and 20 by calculating the difference 
in A^Hi inside and outside the SGSs. The average A^hi in 
SGS 19 and 20 is -3.54 x 10^^ cm-^ and - 3.50 x 10^^ 
cm~^ respectively, while for the surrounding region it 
is -4.90 X 10^^ cm-^ Using the area of SGS 19 and 20 



from jKim et al. ( 1999 ) (circles with diameters of 52.0' and 
49.4' respectively) we estimate the mass lost from SGS 
19 to be -1.6 X 10^ M© and -1.5 x 10^ M© lost from SGS 
20. These masses are remarkably close to the mass of the 
outflow blob. However, this estimate ignores the possi- 
bility that a sizeable fraction of this lost gas has probably 
been swept up by SGS winds and gone into compressing 
the surrounding ISM. Nonetheless, even if only half of 
the lost gas has been blown out of SGS 19 and 20 it is 
sufficient to explain the outflow blob. 

There is another outflow blob of gas to higher ve- 
locity from SGS 1 (Fig. |24|, which is a SGS in arm 
B (the LMC filament) located at (a, ^)=(04^58"^36^ 
-73°33'57'0 with a mean velocity of li sR = 242 km s~\ 
size of (48.6' X 45.4') and age of3.0Myr ("Kim et al.'1999). 
SGS 1 and the outflow can be seen in Figure [24j The 
outflow first shoots to higher velocity by —27 km s~^ at 
nearly the same position and then moves south —0.5°, 
where the densest portion is positioned (Fig. 24 3). The 



blob appears to have two clumps that are separated in a 
by -0.5° (Fig.|24]D) (Fig. [24^). The outflow must have 
substantial tangential motion in order to have moved 
—0.5° (co rresponding to — 0.44 kpc) in —3 Myr, or it is 
older than|Kim et al. ( 1999 ) determined. The mass of the 



< 



blob (within 04^49"^30' <a< 05^05"^48% -74°20'04' 
8 < -73°31'10", 267.4 km s-^ < Flsr < 288.0 km s-^) 
is —7.1 X 10^ Mq, which is about half the mass of the out- 
flow blob coming from SGS 19 and 20. These examples 
of SGSs blowing out massive amounts of gas are strong 
evidence that the SGSs in the SEHO are capable of blow- 
ing out substantial amounts of gas, and, if chained in a 
series of propagated star formation events, creating the 
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Fig. 23. — Integrated intensity of the LMC H I gas from S03 
showing the outflow of gas from SGS 20, centered at (a, (5, Vlsr)= 
(05^43^09% -71°08'46'^ 270 km s-^). (a) Column density, A^m, 
in units of IQ-*^^ atoms cm~^ (integrated from 266 < Vlsr < 
277 km s-i), (b) Vlsr vs. a (integrated from -71°26'34''< S < 
—70° 59^00^^; the greyscale is in units of K), and (c) Vlsr vs. 6 
(integrated from 05^39™48^ < a < 05^44™ 54^; the greyscale is in 
units of K). The locations of SGS 20 and 18 are shown in all three 
panels. Even though the outflow is spatially more aligned with 
SGS 18 than SGS 20 (in panel (a)) it is clear from panels (b) and 
(c) that this is a projection effect. The outflow of gas is connected 
in position and velocity to the bub ble of gas that is surrounding 
SGS 19 and 20 (see Figs. [2l] and [22j, but has moved farther to the 
south and hence is spatially in front of SGS 18. This flgure shows 
that SGS 20 is blowing out a large amount of gas from the LMC. 
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Fig. 24. — Integrated intensity of the LMC H I gas from S03 
showing the outflow of gas from SGS 1, centered at (o;, 5, Vlsr) 
= (04^58^36% -73°33'5r^ 275 km s-^). (a) Column density, 
A^Hl, in units of lO-*^^ atoms cm~^ (integrated along Vlsr for 271 
< ^LSR < 281 km s"-*^), {h) Vlsr vs. ol (integrated along 5 for 
-74°3r33''< 5 < -73°19'59''; greyscale in units of K), and (c) 
^LSR vs. d (integrated along a for 04^57^17^ < a < 05^07^31^; 
greyscale in units of K). This shows that SGS 1 is blowing a sub- 
stantial amount of gas out of the LMC fllament. This also suggests 
that there may be recent star formation going on in the Magellanic 
Stream and there might be some young stars there. 



Leading Arm and the trailing LMC filament of the MS. 

SGS 1 lies wen outside (by -2°) the LMC H I disk 
as defined in Figure [20| and matches arm B perfectl y in 
position and velocity(Fig. [24|. We have argued in ^5.3 
that arm B is the beginning of the LMC filament of the 
MS. If SGS 1 is indeed a part of arm B, as all evidence 
points, and our explanation of arm B is correct, then this 
implies that there is ongoing star formation at the head 
of the MS and that there are probably some young stars 
in the Stream. 

The H I mass of the entire LMC is -4.8x10^ M© (SOS), 
while ~1.3 xlO^ M© of that is in the high- velocity gas. 
The mass of arm B is ~4.9 xlO^ M©, that of arm E 
is -6.9 xlO^ Mo, and -1.2 xlO^ M© is to the west 
of arm B. This high- velocity gas constitutes 27% of the 
measured H I mass of the LMC and —1/3 the combined 
mass of the MS and LAF (which have -3.8 xlO^ M© and 
-3 xlO^ M0 respectively; B03). Thus the LMC has just 
recently blown out and lost a substantial amount of mass 
from the SEHO. It therefore seems plausible that over an 
even more extended period of time the SEHO might have 
produced the Magellanic Stream and LAF. 

What is the mass ledger like for this hypothesis? The 
total mass within the volume of a SGS in the SEHO of 
the LMC with an average diameter of 50' (or 0.72 kpc) is 
- 1 X 10^ M©. If 10% of the mass within a SGS is blown 
out of the LMC disk, then each SGS contributes — 1 x 10^ 
Mq of gas to the MS and/or LAF. If the Magehanic 
Stream originated in the SEHO and was blown out by 
SGSs, then some —410 shells would have been required to 
have created the Stream. If the Stream is —1.74 Gyr old 



(from ^5.5), then the average SGS creation rate would 
have to be 1 per —4.2 Myr. If the shells are coherent 
and observable for —10 Myr, which appears to be th e 
maximum age of shells in the LMC CKim et al | |T999l ), 
then approximately 2-3 should be seen in the SEHO at 
any time. This is the number currently observed (SGS 
18, 19, and 20). 

Of course, one issue with the MS and LAF originating 
in the SEHO is that the current H I mass of the SEHO 
(—1.3x10^ ^0) is now smaller than the H I mass of the 
entire Stream. There are several possible explanations 
for this discrepancy: (1) the SEHO was more massive in 
the past, (2) the SEHO is continually being replenished 
by gas from the LMC disk, or (3) the SEHO is continually 
growing by accreting from the hot, coronal MW it is 
plowing through. Since the current H I masses of the 
LMC and SMC are roughly the same, moving the origin 
of the MS and LAF from the LMC to the SMC does not 
solve the mass problem. The system from which the MS 
and LAF originate has lost a substantial amount of gas 
recently and must have contained at least double the H I 
in the past. 

On the other hand, is it possible for the SGS creation 
rate and g as outflow to be su stained for —1.74 Gyr? Ac- 
cording to |Kim et al.| (|1999), the continual creation of 
SGSs can be acEieved by self-propagating star formation, 
which is observed in the 30 Doradus complex. Therefore, 
it appears plausible for there to have been continuous 
star formation and SGSs blowing out gas sufficient to 
form the MS+LAF for an extended amount of tim e in 
the SEHO. This is discussed further in section \7\2\ 



6.5. Energetics of the Supergiant Shell Blowout 
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Is it energetically realistic to suppose that SGSs have 
blown out the MS and LAF? The energy that would be 
required to blow out the entire mass of the MS and 
LAF (-4.1x10^ M0) to a velocity of -50 km s"^ is 
E = (1/2)Mms^^ ^ 1.0 X 10^^ ergs. Since the average 
kinetic ene rgy released by a SN e is —10^^ ergs ( Woosley| 
& Weaver||1986j^|Cho & Kang||2007|) this would require 



^10,000 SNi: 



Cho & Kang 2007) this would require 
half of the SMe kinetic energy in an 



SGS goes into sweeping up gas and compressing the ISM 
in the disk and the other half blows out gas perpendic- 
ular to the disk, then the combined energy of —20,000 
SNe are required. Spread out over —1.74 Gyr gives an 
average SNe rate of 1 every —87,000 years. This is fairly 
low compared to the MW SNe rate of 1 e very —100 years 
and the LMC rate o f 1 every —500 years (van den BerghJ 
fc Ta m maimlllQQll) . If the SNe are divided into -41(1 
SGSs then each SGS would contain —49 SNe. The 30 
Doradus nebula contains 39 O stars, 12 Wolf-Rayet stars 
and 8 B supe rgiants that will all eventually go supernova 
|Melnick 1985) a tally that validates the plausibility of 
our estimates. Therefore, the energetics appear to be 
favorable to the SGS blowout hypothesis. 

6.6. Magnetic Fields in the SE HI Overdensity 



[Hay nes et al. ( |1991 ) performed a radio continuum sur- 
3y of the Magellanic Clouds at 2.45, 4.75, and 8.55 GHz. 



vey 

In the linearly polarized maps of the LMC Haynes et 
al. discovered that two thin, long (—3-4 kpc) filaments 
(stretching south of the 30 Doradus nebula) dominated 
the emission (most obviou s at 2.45 GHz or 12 cm; see Fig. 



26 1 here or their Fig. 7). |Xu et al.1 ( 19921 ) separated the 
tHermal (free-free) from the nonthermal (synchrotron) 
emission and found that, while there is relatively uni- 
form, nonthermal radiation present across most of the 
LMC (indicating pervasive magnetic fields), the region 
extending south of the 30 Doradus nebula (where the 
filaments are located) contains almost entirely thermal 

emission (see their Figs. 10a and b). 

With t he addition of observations at 1.4 GHz,[Kleinet 



al. ( 1993) used the radio continuum data to study the 3D 
magnetic field structure of the LMC, and found that the 
magnetic field in the region of the filaments was signif- 
icantly aligned (i.e. not random), and oriented more or 
less parallel to the filaments' direction. In addition, the 
most extreme rotation measures (RM) — an indication 
of the magnetic field parallel to the line-of-sight — were 
in this region. Klein et al. conclude that there is a gi- 
ant magnetic loop bending out of the plane of the LMC 
towards us (see their Fig. 6). We compare the high- 
velocity gas with the linearly polarized filaments from 
Klein et al. in Fig. [26^. The two groups of filaments 
(high- velocity H I and "polarized") are offset from each 
other. In fact, they appear to be anti-correlated with 
a low correlation coefficient of 0.30. One exception is 
a small region near SGS 18 (and the SGS 20 outflow 
blob) that is bright in the polarized emission (and in CO 
as weh: Cohen et al. 1988, Fukui et al. 1999), proba- 
bly indicating that this region is very active. It appears 
that the eastern polarized filament might originate from 
there. The western polarized filament might be related 
to the 30 Doradus nebula. The proximity to each other, 
filamentary structure, and abnormality (compared to the 
rest of the LMC) of the two sets of filaments indicates a 
probable physical relationship between them; the exact 
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Fig. 25. — Velocity slices of the H I datacube from S03 showing 
where the filaments of arm B (the LMC filament) and arm E (the 
LAF) originate (the greysc ales are in units of K). The outlines of 
SGS 1, 18, 19 and 20 from |Kim et al. |(|1999| are also shown. In 
panel (a) at Vlsr — 190.7 km s~^, the two southern filaments of 
arm B (with fairly faint Tb) appear to originate in SGS 18. In 
panel (6) at Vlsr = 255.0 km s"-*^, the northern filament of arm 
B and a part of arm E seem to originate in SGS 18. Panel (c) at 
Vlsr — 273.2 km s"-*^, the northern filament of arm B emanates 
from the north-western part of SGS 18. Dense knots of H I gas 
around SGS 19 and 20 seem to be the source of some of the arm 
E gas. 
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Fig. 26. — The sky distribution (c olum n density, A^hi? in units of 10-*^^ atoms cm"-^) of the high- velocity component of the LMC H I gas 
from the S03 datacube (same as Fig.[20p). (a) The 12 cm linearly polarized radio emission from Klein et al. ( 1993) is overplotted as cont our 
lines and shows the two magnetic "tilaments" . There is strong emission from 30 Doradus, the ISGIS 20 outflow blow (also see Fig. |23| >, as 
well as two filaments which are parallel to the high-velocity gas, especially arm B. The eastern polarized filament appears to be associated 
with the SGS 20 outflow blob. The polarized and H I sets of filaments seem to be anti-correlated and a small correlation coefficient of 0.30 
bears this out. (b) The density of 2MASS LMC red giant branch stars (Skrutskie et al. 2006) is overplotted with contour lines showing the 
central stellar LMC bar. The eastern end of the bar is near the region where most of the high- velocity gas is originating. The proximity of 
the end of the bar might be partly responsible for the high-density of H I gas and perturbed dynamics in this region of the LMC. 



nature of that relationship and the origin of the spatial 
offset remain unclear at this time. 

These polarization and magnetic-field studies show 
that there is something violent happening in the region 
south of 30 Doradus that is likely related to the high- 
velocity gas and the SGSs. This corroborates circum- 
stantially our suggestion that the SGSs may be blowing 
out the H I gas and creating the LMC filament of the 
Magellanic Stream and the Leading Arm. It is possi- 
ble that the relativistic ions are being accelerated by the 
supernovae shocks from the SGSs, and that the aligned 
magnetic field is being caused by the stellar winds and 
supernovae shocks (also from the SGSs). It is not en- 
tirely clear how the magnetic loop fits into this picture. 
However, if there are magnetic fields in the Magellanic 
Stream it might explain how the filaments can remain so 
coherent over large distances. 

6.7. Relative Distance of the High-Velocity Gas 

There is some diagreement about whether the anoma- 
lous ( ''hig h- velocity" ) ^as is in front of or behind the 



tures that correlated spatially with the distribution of 
the high-velocity gas. Based on this evidence they con- 
cluded that the high- velocity gas is behind the plasma in 
LMC 2 (which presumably is in the disk of the LMC). 
However, this study only looked at a small subsection of 

t he LMC. 

Dickey et al.|(1994) performed an absorption-line study 



LMC disk. |Luks &; Rohlfe] ( |1992D conclude that the low- 
velocity component is in front of the high-velocity com- 
ponent since they find no absorption of the 30 Doradus 
emission (associated with the low - velocity com ponent) 



by the high- velocity component. Points et al.| ( |1999[ ) 
used ROSAT data to study the supergiant shell LMC 
2 (SGS 19). They detected no X-ray absorption fea- 



of several sight-lines in the LMC. They used a kinemat- 
ical argument to conclude that the high- velocity compo- 
nent is in front of the low- velocity component. Since 
the LMC disk is inclined, a cloud at a height above 
the LMC disk plane (and following the same rotation 
curve) will have an actual LMC-centric distance that is 
different than the projected one and therefore will have 
a different rotational velocity than the rest of the gas 
along that same sight-line (that is in the disk). Whether 
the rotation velocity will increase or decrease depends 
on the exact sight-line, and this information is embodied 
in a "velocity gradient" ((i'i;/(ir). Dickey et al. calculated 
the velocity gradient for each position at which they de- 
tected absorption lines using the smoothed velocity field 
by Luks & Rohlfs (1992; their Fig. 9) and assuming that 
the north-east is the near side of the LMC. Based on the 
velocity gradient at each position and the radial velocity 
offset from the disk radial velocity at that point, they 
conclude that the low-velocity gas is on the far side of 
the disk and the high- velocity gas is in front. 
The kinematical argument by Dickey et al. works fairly 
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well for arm B, which moves across a large part of the 
LMC. Based on the velocity gradient its velocity offset 
should at first be positive, then zero, and then negative, 
and this is actually what is observed. However, there is 
a problem with the gas "blobs" that are being blown out 
of the SGS mentioned above (^6.2). The gas moves to 
higher velocity at nearly the same position which would 
be interpreted to mean, without considering the rotation 
curve, that the gas is moving away from us. Furthermore, 
since it seems clear now that most, if not all, of the high- 
velocity gas in the south-eastern portion of the LMC is 
coming from the SGSs, the high-velocity gas should be 
on the same side of the disk as the "blob" gas. Even 
though the kinematical argument of Dickey et al. seems 
realistic, it assumes that the gas out of the disk plane 
is also in circular motion. However, we question if this 
assumption is a valid one. It is clear that the LMC's 
gravity is not the only force acting on the gas. As al- 
ready stated, it seems clear that the SGSs are blowing 
gas out of the disk (probably away from us), and since 
arm B becomes the ~100° trailing LMC filament of the 
MS, tidal or ram pressure forces must be acting on it. 
Therefore, it seems unlikely that a simple model with the 
gas in relative equilibrium, such as proposed by Dickey et 
al., will work correctly in this highly dynamical situation 
with various contributing forces. 



ig lorces. 

( [l993| ) claim that two filaments 
diarized emission are on the near- 



Finally, Klein et al. 
visible in the linearly poli 
side of the disk. They conclude that the filaments are 
outside the disk since their rotation measures are much 
lower than would be expected if they were inside the 
LMC disk. The rotation measures would also be higher 
than expected if the filaments were on the far side of the 
disk. Due to the anti-correlation of the high- velocity H I 
and "polarized" filaments (mentioned above) it is doubt- 
ful that the claims by Klein et al. about the polarized 
filaments can be used to ascertain the position of the 
high-velocity H I gas relative to the disk. However, the 
possible association of the H I filaments with the polar- 
ized filaments, and the "magnetized loop" hypothesized 
by Klein et al., deserve further study. 

We conclude that the high-velocity gas of both arms 
E and B is on the far side of the LMC, mainly because 
we are attributing the dynamics to the force of the SGSs 
blowing out gas to higher velocity, and this must put 
higher- velocity gas away from us. This is consistent with 
the conclusions of both Luks & Rohlfs and Points et al., 
but inconsistent with Dickey et al. as well as with the 
Klein et al. assessment of the polarized filaments. 

7. DISCUSSION AND SUMMARY 

Our exploration of the Magellanic System using a 
Gaussian decomposition of the LAB data in combina- 
tion with other radio data has led us to several conclu- 
sions regarding the nature and origin of the Magellanic 
Stream. 

7.1. The Large Magellanic Cloud as Progenitor of the 
Magellanic Stream and Leading Arm Feature 

We have found evidence that one of the filaments of the 
Magellanic Stream and the Leading Arm can be traced 
back to the LMC and that both of these H I structures 
originate there. In our database of H I Gaussians the 
MS appears as two filaments as has been previously re- 



ported by [Putman et al." (2003). But, capitalizing on 
the coherence of the filaments m ppv (position-position- 
velocity) space, we were able to track one of them back 
to the LMC in both velocity and position. After isolat- 
ing the H I of this filament in one of the position- velocity 
diagrams (Flsr vs. Lms)^ we were able to show that it 
originates in the SEHO of the LMC. Likewise we showed 
that the first complex of the LAF (LA I) begins (in posi- 
tion and velocity) near the end of S03's arm E from the 
LMC. The spatial and velocity progression of the three 
concentrations of LA I and arm E, as well as their conti- 
nuity in the deeper HIPASS data, strongly suggest that 
they are physically connected. 

An LMC origin of the MS is contrary to much of 
the current literature, which has largely supported an 
SMC/Bridge origin for the MS (e.g., P03 and BOS). Even 
most modeling papers (especially those based on a tidal 
origin of the MS) have started with the SMC/Bridge as- 
sumption and adopt an N-body representation of the 
SMC influenced by only a static LMC (and MW) po- 
tential (e.g., Gardiner & Noguchi 1996; Yoshizawa & 
Noguchi 2003; Connors et al. 2004, 2006). Therefore, 
these models by design rely on the MS forming from the 
SMC and not the LMC. Models that represent both the 
LMC and SMC as N-bodies (Murai & Fujimoto 1980; 
Ruzicka et al. 2006) have also concluded that most of 
the Magellanic Stream came from the SMC; however, 
these simulations do show at least some debris coming 
off the LMC and this debris follows the position and ve- 
locity distribution of the main body of the Stream. The 
main reason that an SMC origin was preferred is because 
more material came off the SMC than the LMC in these 
models. The reasons cited for this observed mass-loss 
inequity is that the SMC potential is weaker than the 
LMC's and that the LMC has a much larger tidal influ- 
ence on the SMC than vice versa. Even though these 
statements may be correct in the context of these mod- 
els, it is likely that these models are missing important 
physical processes: none of the simulations include the 
potentially important dynamics of LMC star formation 
and SGSs that may greatly enhance gas outflow from the 
LMC. If SGSs propel gas faster than the escape velocity, 
this gas can then be swept away from the MCs by tidal 
and/or ram pressure forces. Future simulations includ- 
ing these effects subsequent to blowout may come closer 
to matching the H I structural features that we have ob- 
served. First steps towards this goal have already been 
taken by Olano (2004) who has already demonstrated 
that blowout from the MCs can create the large scale 
features of the MS, as well as the HVC system of the 
MW. Olano's model explored blowout from the MCs that 
results from an interaction between the two Clouds 570 
Myr ago. Our analysis here identifies a more specific 
site (the SEHO of the LMC) and timescale (-1.7 Gyr) 
for the creation of the MS that we believe ought to be 
accommodated by such models. 

With the discovery of the two filaments of the MS 
(P03) it was suggested that one of the filaments orig- 
inates from the Bridge, and that therefore the Bridge 
migh t be older than the —200 Myr s uggested by the mod- 
els of Gardiner, Sawa, fc Fujimotol fl994) and the young 
blue stars discovered between the Magellanic Clouds by 



Irwin, Demers fc Kunkel^ (1990) . However, now that this 
same H I filament can be traced to the LMC there is 
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no reason to believe that the Bridge is old. In fact , the 
rece nt stu dy of stellar populations in the Bridge by |Har 



(2006) shows that there are only young stars there 
and that star formation started in the Bridge only some 
200-300 Myr ago. 

According to Russell & Dopita (1992), the current 
mean metallicities of gas in the LMC and SMC are 
[Fe/H]= -0.2±0.2 and -0.6±0.2, respectively, while 
the average metallicity of the MS is [Fe/H]= -0.6±0.2 
(Wakker 2001). At first glance it appears that the metal- 
licity information points towards an SMC origin of the 
MS. However, it is not the current metallicity of the MCs 
that is important, but the metallicity at the time that 
each part of the MS left the Clouds, starting around 
1.74 Gyr ago (so that one might expect a metallic- 
ity gradient along the Stream) . According to ^Pagel fc 
|Tautvaisiene| ( 1998 ) , the metallicities of the Clouds were 
^^0.3 dex lower ~1.7 Gyr ago, which would put them at 
[Fe/R]^ -0.5 (LMC) and [Fe/R]^ -0.9 (SMC). Thus, 
even though metallicity provides only a weak discrimi- 
nant of the MS origins, it appears to slightly favor an 
LMC origin over an SMC origin. 

Moreover, there continues to be a global problem rec- 
onciling large observed variations in the metallicities of 
different gaseous Magellanic features. For exam ple, the 
low metallicity of the Bridge, [Fe/H]= —1.1 (Lehner 
|2002), compared to [Fe/H]= -0.6 for the SMC (where 
it is thought to have originated), has been used to argue 
that the Bridge is old. However, because other evidence 
supports a young age (~200 Myr) for the Bridge (see 
above), the enrichment level of the gas must be more 
complex than mere mass loss fr om a simple c losed box 
model of the MCs. For example, Lehner (2002) suggests 
that the Bridge gas from the SMC may have mixed with 
other, less enriched gas. One possible source of less en- 
riched gas could be earlier gas ejections from the SMC 
itself which could be accreted. Another possible expla- 
nation for the low metallicity of the Bridge gas is that 
it could have originated from gas in the outskirts of the 
SMC that might be more metal-poor than gas in the 
central part of the SMC. Whatever the explanation is, it 
seems clear that it is difficult at present to interpret the 
relative metallicities of Magellanic gas and the MCs. 

7.2. Origins in the SE HI Overdensity 

We have been able to track the MS back to the LMC, 
and argue, in ^ for its origin in the SEHO. What re- 
mains unresolved is whether the physics creating this in- 
tense star- forming region is tied to a particular, fixed 
dynamical hot spot that is the product of the bulk mo- 
tion of the LMC, its rotation, and interaction with MW 
halo gas, or whether the SEHO participates in the gen- 
eral LMC rotation. The answer bears on how to evaluate 
an accounting of relative ages, masses, and energetics be- 
tween the SEHO site and what is needed to produce the 
MS and LAF an d their apparent oscil lations. 

For example, IdeBoer et al. fl998") suggest that the 
leading edge of tIie~LMC (the SE, where the SEHO 
is presently located) is being compressed as the LMC 
moves through the diffuse halo of the MW (creating a 
bow-shock), which produces a steep H I gradient (as ob- 
served), and triggers star formation. In this scenario, the 
southeast quadrant of the LMC will always remain an ac- 
tive star- forming site, with previous generations of stars 



formed there rotating off of the hot-spot due to LMC ro- 
tation, but with the evolution of massive stars constantly 
regenerating SGSs near their birthplace. This argues for 
considering a broader perspective in evaluating the age 
and energetics of MS production. Unfortunately, while 
it is obvious that there has been much recent star forma- 
tion occuring in the 30 Doradus (Tarantula) nebula and 
its immediate surroundings, the star formation history of 
the SEHO and the annulus of the LMC it occupies have 
been less well established. But from our assessment of 
the mass, age, an d req uire d en ergy needed to produce the 
MS and LAF (§ ^6.4| and 6.5) we can hypothesize some 
general expectations. 

We have given evidence in ^that SGSs in the SEHO 
are currently blowing out gas and creating the MS and 
LAF. This scenario is supported by the high-speed out- 
flows from the 30 Doradus starburst ( j Redman et al.|20"Q3]) 
and by the energetic outflows from the LMC ( [Lehner 
Howk|20'Q7 ). From the continuity of the MS we can infer 
that tnis process must have been going on continuously 
at a rate steady enough to not create gaps in the MS. 
In the de Boer et al. scenario one might expect a contin- 
uous cycle of star formation, SGSs, and gas blowout as 
long as the compression and bow-shock are present. As 
discussed in §6.3, to eject the -4.1x10^ M© of the MS 
and LAF out of the SEHO in ~1.7 Gyr would require 
-400 SGSs (each blowing out -10^ M©). This exceeds 
the current H I mass of the SEHO (-8x10^ M©, exclud- 
ing the high- velocity gas), but this comparison belies a 
complex calculus of mass exchange. On the one hand, 
some H I from the SEHO is converted into stars and 
stellar remnants, while some is lost to SGS blowout to 
produce the MS and LAF, and additional H I gas can 
be generated from the destruction of H2, although this 
is likely to be a small amount.^ On the other hand some 
gas may be inflowing from the MW halo (and indeed may 
be the source of less enriched material complicating the 
metallicity accounting discussed above). 

If the MS and LAF were formed as part of the de Boer 
et al. scenario, then a >1.7 Gyr age for these gaseous 
structures is viable, but the star formation that pro- 
duced these features would have created stellar popu- 
lations that would have rotated around the LMC about 
five times. Therefore, we would expect that the stellar 
populations at any position in a 2 kpc radius to have five 
distinct peaks in age corresponding to the times when 
this position was in the SE corner of the LMC where 
the vigorous star formation is occurring. These peaks 
should be separated by —340 Myr (i.e. the LMC rota- 
tion period at this radius) but with approximately the 
same star formation rate. The exact ages of these peaks 
should be shifted with position in the 2 kpc annulus - 
older as you move clockwise (with the LMC rotation) 
and younger in the opposite direction. The total number 
of stars younger than —1.7 Gyr should be approximately 
constant with position around the 2 kpc annulus. Thus, 
a detailed analysis of the star formation history of the 2 
kpc annulus could be used to check the validity of the de 
Boer et al. scenario. Unfortunately, radial and azimuthal 
mixing of stars in the LMC disk might smear out some 



^ According to [Young ScoviHe| ( |l99l| the mass in molecular 
hydrogen is typically an order ot magnitude less than that in atomic 
form for very late-type spirals. 
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of these patterns which could make this analysis difficult 
in practice. 

It is not clear what caused the vigorous star formation 
in the SEHO region to start ~1.7 Gyr ago and com- 
mence the SGS blowout and formation of the MS and 
LAF. One possibility is that the LMC and SMC had 
a close encounter that could have triggered star forma- 
tion. The inte grated star formation histories of the LMC 
(iHarris & Zaritsky 2004) and SMC ( Smecker-Hane et al. 
^ both show 



a "burst" of star formation ~2 2.5 Gyr 
ago which suggests that an interaction between the MCs 
might indeed have occurred around that time. The LMC 
globular cluster age-gap between 3 and 13 Gyr (e .g., Da 
T99TI IGeisler et aL][T997| [Rich et al.||2QQl| |Piatti 



Costa 



et al.| 2002 ) indicates an onset of star formation in the 
LMC around ~3 Gyr ago, h owever, such an age-gap is 
not seen in the SMC clu sters ( Da Costa|1991|). As prevr 



ously mentioned in ^5.5, the orbits by Besla et"aL] (2007) 
indicate possible encounters of the MCs with eacn other 
around 0.2, 3 and 6 Gyr ago. Therefore, it does not seem 
unrealistic to suppose that a close encounter of the MCs 
~2-3 Gyr ago triggered star formation in the LMC and 
SEHO region. Another possibility might be that ~1.7 
Gyr ago the LMC was close enough to the MW that ram 
pressure due to the hot MW halo created a bow shock 
and started star formation in the SEHO (jump starting 
the de Boer et al. scenario mentione d above). However , 
according to the hyperbolic orbits of |Besla et al.| p007 ) 
the LMC was -300-500 kpc from the MW at This time 
(depending on the MW mass); even though the density 
of the gaseous MW halo is not well known it is doubtful 
that the density at those distances could be large enough 
to create the bow shock necessary to start star formation 
in the LMC. The trigger of the recent, vigorous star for- 
mation in the SEHO region of the LMC therefore remains 
an open, and interesting question. 

It should be acknowledged that the SMC also has many 
giant shells and SGSs (Staveley-Smith etaL 1997), and it 
might be expected that the SMC is also blowing out gas 
that might be contributing to the MS. For example, it 
may be that the second filament of the MS is coming from 
the SMC. Unfortunately, we have no evidence to support 
or disprove this based on our analysis. The lack of an 
obvious major source of frenetic star formation activity 
in the SMC may substantially lower the efficiency of and 
amount of SMC blowout relative to that in the LMC. 
Whether the SMC contributes to the MS and by how 
much requires further work to determine. 

7.3. The Sinusoidal Velocity Pattern of the LMC 
Filament 

Among the most striking features of the MS, as seen in 
Figures [9]3 and[Toj are the velocity oscillations of the two 
filaments. These sinusoidal patterns have not been previ- 
ously observed or predicted. We have hypothesized that 
this pattern for the LMC filaments may be an imprint of 
the motion of the gas ejection site according to the LMC 
rotation curve. Using this hypothesis, we estimated that 
the drift rate of the MS gas away from the LMC is —49 
km s~^, and, based on the length of the entire MS, we 
surmise that the MS is —1.7 Gyr old. As we mentioned 
in the previous subsection, this estimate of the age of the 
MS is fairly consistent with other studies. 

However this hypothesis contradicts the de Boer et 



al. (1998) scenario by invoking the SEHO as a rotat- 
ing star- formation site. In this case, the multiple gen- 
erations of stars formed in this process should still be 
relatively near the ejection site. But without replenish- 
ment of gas, this scenario may require inordinately large 
amounts of mass to be processed through one particu- 
lar star-formation site with no apparent driver and for a 
relatively long time. 

A hybrid hypothesis that could incorporate the driving 
physics of the de Boer et al. scenario and still produce 
oscillating patterns in the MS and LAF might include gas 
streaming along the LMC bar, creating a "perfect storm" 
of compressional activity w hen combined with bulk LMC 
motion and LMC rotation. |Kim et al. ( 1998 ) show that 
the LMC H I velocity field deviates from circular rota- 
tion, especially at the north-western end of the stellar 
bar, and take this as evidence for large-scale streaming 
motions in the H I. In Figure [26) 3 we show the high- 
velocity H I gas together with the density of 2MASS 
LMC red giant branch stars (Skrutskie et al. ^06) indi- 
cating the stellar bar. The eastern end of the bar is close 
to the SEHO where the H I filaments originate, suggest- 
ing that the bar may play a role in (1) the accumulation 
of the dense gas in SEHO and (2) the origin of the high- 
velocity filaments there. However, whereas the end of 
the bar on the leading side of the LMC coincides with in- 
tense star formation in the SEHO, no equivalent feature 
is found on the western side of the bar. If bar streaming 
provides the critical, extra dynamical contribution nec- 
essary for star formation activity at the level needed for 
persistent blowout, then oscillating motion of the ejec- 
tion site would naturally correspond to the bar pattern 
rotation (which is not well known). But this would also 
imply periodic jumps in location every half period when 
there is a transition from one end of the bar being on the 
leading side of the LMC to the other end leading. This 
motion would be imprinted in the shape and velocity of 
the MS. Spatial jumps similar to that suggested may be 
visible (e.g.. Fig. 13), but they are not as evident in the 
velocity distribution. 

Clearly hydrodynamical modeling is needed to resolve 
the dilemma of how the ejection site is generated and 
to understand the expected effects of its motion, or lack 
thereof, on the MS and LAF. We consider this one of the 
most important challenges remaining to complete this 
new picture of how the MS and LAF formed. 

7.4. Relevance to the Tidal vs. Ram Pressure Models 

The MS origins debate has focused on the tidal versus 
ram pressure models. We have added a new mechanism, 
SGS blowout, to explain how the MS gas is removed from 
the MCs. Once the gas has escaped its host galaxy, ram 
pressure and/or tidal forces are still required to disperse 
the unbound gas and move part of it forward (to cre- 
ate the LAF) and other parts backward (to create the 
MS). Therefore, the blowout model moves the tidal ver- 
sus ram pressure debate from the mechanism for removal 
of gas from the MCs, to the mechanism for dispersal of 
the unbound gas. 

It seems likely that both ram pressure and tidal forces 
are present and needed to explain the characteristics of 
the MS and LAF: The LAF cannot be satisfactorily ex- 
plained without a tidal force, and ram pressure forces are 
evidently at work in building the steep density gradient 
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in the leading edge of the LMC. The combination of ram 
pressure and tidal forces can also help explain t he column 
density ffladient alon^ the Magellanic Stream (Moore & 
Davis||l994 Mastropietro et al. 2005 ) and the imbalance 



of mass in the MS and LAF, since it will be more diffi- 
cult to move material ahead of the MCs due to the extra 
force pushing it backwards. 

Our proposed scheme for the origin of the MS through 
blowout resolves a problem that has plagued tidal mod- 
els, namely the lack of observed stars in the Stream. If 
most of the MS and LAF gas was blown out of the LMC 
from SGSs in the SEHO, as suggested here, then no stars 
would be expected in the Stream. The forces operating 
in the SGSs that blow out the gas (superwind and super- 
novae shocks) do not affect the stars, and therefore the 
mystery of the lack of stars in the Stream, even in the 
presence of tidal mechanisms, is easily explained. 

Our mode l also r esolves the paradox recently posed by 
[Besla et aT] ( |2007| ) that if the MCs are on hyperbolic 
orbits, as now in dicated by HST proper motions (Kalli- 
l^ayalil et al. 2 006a|b} [Piatek, Pryor, fc Olszewski 2007), 
then the distance of the MCs from the MW at the time 
that the MS originated is too large (^300-500 kpc) for 
the ram pressure and tidal forces to strip the gas out 
of the LMC or the SMC. SGS blowout can do the work 
of moving the gas to large enough LMC radii to where 
even weak ram pressure and tidal forces can take over for 
dispersal in the hyperbolic orbit scenario. 

7.5. New Constraints For Modeling of the Magellanic 

Stream 

In this paper we have used the detailed spatial and 
velocity distributions of Magellanic H I gas to lead us 
to a new paradigm for how the MS and LAF formed, 
namely through blowout from intense star formation in 
the SEHO of the LMC. But there are intriguing patterns 
in the MS and LAF that can provide additional clues at 
a more detailed level to the processes that shaped them. 
The most obvious clues are the periodic patterns. Among 
them: 

1. The two MS filaments exhibit large velocity and 
spatial oscillations. The oscillation of the LMC fil- 
ament has a velocity amplitude of 26.4 km s~^ with 
an angular period of 19.5° (17.1 kpc at a distance 
of 50 kpc) in Lms (Fig. [17|). The amplitude of the 
spatial oscillation in Bus is ~2° (Fig. 13). 



2. There are three concentrations in LA I, elongated 
along Bus ("^2 x 7°) that look very similar to 
one another and that are each offset by ~12° in 



Lms and ~9.5° in Bus from the previous one (Fig. 



3. LA II and LA III are very similar in appearance 
(Fig. Is]). They are both elongated along Lms with 
sizes of ~21 X 5° and parallel to each other (offset 
by -20° in Bus)- 

4. The two filaments of the Magellanic Stream exhibit 
strong periodic patterns in position, and are com- 
posed of clumps elongated along Lms with sizes of 
—6.0x1.5°, surprisingly similar in size to the LA I 
clumps and highly suggestive that they may have 
formed by the same process (Fig.js})). 

5. The two MS filaments are quite similar in appear- 
ance and mirror each other in their shape (from 
Lms - -15° to -45°), only shifted by -1° in 
Lms and -4° in Bus (Figs.jSjD and[l3). 



Future modeling efforts of the MS should not only incor- 
porate an LMC/SEHO origin for both the MS and LAF, 
but strive to reproduce these other newly found distinc- 
tive observational characteristics of the MS. We propose 
that such models account for the energetics of supergiant 
shell blowout in the creation of the MS. 



We appreciate useful discussions with Bill Kunkel, 
Gurtina Besla, Pavel Kroupa, Jon Hibbard, Remy Inde- 
betouw, Ed Murphy, Zhi-Yun Li, Jake Simon, John Haw- 
ley and Amy Reines. Our work has made use of the IDL 
program MPFIT by Craig Markwardt, and we would like 
to thank him for making this and other useful programs 
available to the public. We would also like to thank Lis- 
ter Staveley- Smith for providing us with the Parkes H I 
datacube of the LMC as well as Bryan Gaensler and Uli 
Klein for providing us with the 12 cm polarization data of 
the LMC. D.L.N, is supprted by the ARCS Foundation, 
the Green Bank Telescope Student Support Program, 
a University of Virginia President's Fellowship, and the 
Virginia Space Grant Consortium. S.R.M. acknowledges 
funding from NSF grant AST-0307851, NASA/JPL con- 
tract 1228235, and the generous support of Frank Levin- 
son and through the Celerity Foundation. This publica- 
tion makes use of data products from the Two Micron 
All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analy- 
sis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and 
the National Science Foundation. We would like to thank 
the referee Jacco van Loon for numerous useful comments 
and suggestions. 



REFERENCES 



Arnal, E. M., Bajaja, E., Larrarte, J. J., Morras, R., &; Poppel, W. 

G. L. 2000, A&AS, 142, 35 
Bajaja, E., Cappa de Nicolau, C. E., Martin, M. C, Morras, R., 

Olano, C. A., & Poppel, W. G. L. 1989, A&AS, 78, 345 
Bajaja, E., Arnal, E. M., Larrarte, J. J., Morras, R., Poppel, W. 

G. L., &; Kalberla, P. M. W. 2005, A&A, 440, 767 
Barnes, D. G., et al. 2001, MNRAS, 322, 486 

Beichman C., Neugebauer G., Habing H. J., Clegg P. E., Chester 
T. J., eds., 1988, IRAS Catalogues and Atlases Explanatory 
Supplement, NASA RP-1190, Vol. 1. GPO, Washington, DC 

Belokurov, V., et al. 2006, ApJ, 642, L137 

Besla, G., Kallivayalil, N., Hernquist, L., Robertson, B., Cox, T. J., 
van der Marel, R. P., & Alcock, C. 2007, ApJ, 668, 949 



Bica, E.L.D., Schmitt, H.R., Dutra, CM., & Oliveira, H.L. 1999, 

A&A, 117, 238 
Bolatto, A. D., et al. 2007, ApJ, 655, 212 
Briick, M. T., & Hawkins, M. R. S. 1983, A&A, 124, 216 
Briins, C, et al. 2005, A&A, 432, 45B (B05) 
Briins, C, Kerp. J., &; Pagels, A. 2001, A&A, 370, L26 
Bullock, J. S., & Johnston, K. V. 2005, ApJ, 635, 931 
Burton, W.B. 1970, A&AS, 2, 291 

Cho, H., & Kang, H., 2007, New Astronomy, submitted 

(arXiv/0708.3478) 
Cohen, R.J. 1982, MNRAS, 199, 281 

Cohen, R. S., Dame, T. M., Garay, G., Montani, J., Rubio, M., & 
Thaddeus, P. 1988, ApJ, 331, L95 



30 



Origin of the Magellanic Stream 



Connors, T. W., Kawata, D., Maddison, S. T., & Gibson, B. K. 

2004, PASA, 21, 222 

Connors, T. W., Kawata, D., & Gibson, B. K. 2006, MNRAS, 371, 
108 

Da Costa, G. S. 1991, in lAU Symp. 148, The Magellanic Clouds, 
ed. R. Haynes & D. Milne (Dordrecht: Kluwer), 183 

Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 1985, 
ApJ, 292, 371 

de Boer, K. S., Braun, J. M., Vallenari, A., & Mebold, U. 1998, 
A&A, 329, L49 

Dickel, J. R., Mclntyre, V. J., Gruendl, R. A., & Milne, D. K. 2005, 
AJ, 129, 790 

Dickey, J. M., Mebold, U., Marx, M., Amy, S., Haynes, R. F., & 

Wilson, W. 1994, A&A, 289, 357 
Epchtein, N., et al. 1997, Messenger, 87, 27 
Fukui, Y., et al. 1999, PASJ, 51, 745 

Gardiner, L. T., Sawa, T., &; Fujimoto, M. 1994, MNRAS, 266, 567 

Gardiner, L. T., & Noguchi, M. 1996, MNRAS, 278, 191 

Geisler, D., Bica, E., Dottori, H., Claria, J. J., Pitatti, A. E., & 

Santos, J. F. C, Jr. 1997, AJ, 114, 1920 
Grillmair, C. J. 2006a, ApJ, 645, L37 
Grillmair, C. J. 2006b, ApJ, 651, L29 
Grillmair, C. J., & Dionatos, O. 2006, ApJ, 641, L37 
Guhathakurta, P., & Reitzel, D. B. 1998 in ASP Conf. Ser. 136, 

Galactic Halos: A UC Santa Cruz Workshop, ed. D. Zaritsky 

(San Francisco: ASP), 22 
Hartmann, D., & Burton, W. B. 1997, Atlas of Galactic Neutral 

Hydrogen, Cambridge University Press, ISBN 0521471117 
Harris, J. 2007, ApJ, 658, 345 
Harris, J., & Zaritsky, D. 2004, AJ, 127, 1531 
Hand, U. 2000, A&A, 364, 83 
Haynes, R. F. 1991, A&A, 252, 475 
Heller, P., &; Rohlfs, K. 1994, A&A, 291, 743 
Hibbard, J. E., Vacca, W. D., & Yun, M. S. 2000, AJ, 119, 1130 
Ibata, R., Lewis, G. F., Irwin, M., Totten, E., & Quinn, T. 2001, 

ApJ, 551, 294 

Irwin, M. J., Demers, S., & Kunkel, W. E. 1990, AJ, 99, 191 
Johnston, K. V. 1998, ApJ, 495, 297 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., Arnal, E. M., 
Bajaja, E., Morras, R., & Poppel, W. G. L. 2005, A&A, 440, 775 
Kalberla, P. M. W., & Hand, U. 2006, A&A, 455, 481 
Kallivayalil, N., van der Marel, R. P., & Alcock, C. 2006b, ApJ, 
652, 1213 

Kallivayalil, N., van der Marel, R. P., Alcock, C, Axelrod, T., 
Cook, K. H., Drake, A. J., & Geha, M. 2006b, ApJ, 638, 772 

Kaper, H. B., Smits, D. W., Schwarz, U., Takakubo, K., & van 
Woerden, H. 1966, Bull. Astron. Inst. Netherlands, 18, 465 

Kim, S., Dopita, M. A., Staveley-Smith, L., & Bessell, M. S. 1999, 
A&A, 118, 2797 

Kim, S., Staveley-Smith, L., Dopita, M. A., Freeman, K. C, Sault, 
R. J., Kesteven, M. J., &; McConnell, D. 1998, ApJ, 503, 674 

Kim, S., Staveley-Smith, L., Dopita, M. A., Sault, R. J., Freeman, 
K. C, Lee, Y., & Chu, Y.-H. 2003, ApJS, 148, 473 

Klein, U., Haynes, R. F., Wielebinski, R., & Meinert, D. 1993, 
A&A, 271, 402 

Konz, C, Briins, C, & Birk, G. T. 2002, A&A, 391, 713 

Kunkel, W. E., Irwin, M. J., &; Demers, S. 1997, A&AS, 122, 463 

Lehner, R. 2002, ApJ, 578, 126 

Lehner, R., Sz Howk, J. C. 2007, MNRAS, 377, 687 
Lin, D. N. C, &; Lynden-Bell, D. 1977, MNRAS, 181, 59 
Lin, D. N. C, &; Lynden-Bell, D. 1982, MNRAS, 198, 707 
Luks, T., & Rohlfs, K. 1992, A&A, 263, 41 

Majewski, S. R., Skrutski, M. F., Weinberg, M. D., &; Ostheimer, 

J. C. 2003, ApJ, 599, 1082 
Mastropietro, C, Moore, B., Mayer, L., Wadsley, J., & Stadel, J. 

2005, MNRAS, 363, 509 

Mathewson, D. S., Cleary, M. N., & Murray, J. D. 1974, ApJ, 190, 
291 

Mathewson, D. S., Ford, V. L., Schwarz, M. P., & Murray, J. D. 
1979, in Proc. lAU Symp. 84, The Large-Scale Characteristics 
of the Galaxy, ed. W.B. Burton (Dordrecht: Reidel), 547 

McGee, R. X., & Newton, L. M. 1986, PASA, 6, 471 

Meaburn, J. 1984, MNRAS, 211, 521 

Meaburn, J., & Blades, J. C. 1980, MNRAS, 190, 403 

Meaburn, J., Marston, A.P., McGee, R. X., & Newton, L. M. 1987, 
MNRAS, 225, 591 

Meixner, M., et al. 2006, AJ, 132, 2268 

Melnick, J. 1985, A&A, 153, 235 



Meurer, G. R., Bicknell, G. V., & Gingold, R. A. 1985, PASA, 6, 
195 

Mihos, J. C. 2001, ApJ, 550, 94 

Mill, J. D., O'Neil, R. R., Price, S., Romick, G. J, Uy, O. M., & 
Gaposchkin, E. M. 1994, J. of Spacecraft & Rockets, 31, 900 

Mizuno, N., et al. 2001, PASJ, 53, 971 

Moore, B., & Davis, M. 1994, MNRAS, 270, 209 

Moore, B., Ghigna, S., Governato, F., Lake, G., Quinn, T., Stadel, 
J., & Tozzi, P. 1999, ApJ, 524, L19 

Morras, R. 1982, A&A, 115, 249 

Morras, R. 1983, AJ, 88, 1 

Muller, E., Staveley-Smith, L., Zealey, W., & Stanimirovic, S. 2003, 

MNRAS, 339, 105 
Munoz, R. R., et al. 2006, ApJ, 649, 201 
Murai, T., & Fujimoto, M. 1980, PASJ, 32, 581 
Navarro, J. F., Frenk, C. S., &; White S. D. M. 1996, ApJ, 462, 563 
Navarro, J. F., Frenk, C. S., &; White S. D. M. 1997, ApJ, 490, 493 
Newberg, H. J., et al. 2002, ApJ, 569, 245 
Odenkirchen, M., et al. 2001, ApJ, 548, L165 
Olano, C. A. 2004, A&A, 423, 895 

Pagel, B. E. J., & Tautvaisiene, G. 1998, MNRAS, 299, 535 
Philip, A. G. D. 1976a, BAAS, 8, 532 
Philip, A. G. D. 1976b, BAAS, 8, 352 

Piatek, S., Pryor, C, & Olszewski, E. W. 2007, AJ, accepted 
Piatti, A., Sarajedini, A., Geisler, D., Bica, E., & Claria, J. J. 2002, 

MNRAS, 329, 556 
Points, S. D., Chu, Y. H., Kim, S., Smith, R. C, Snowden, S. L., 

Brandner, W., & Gruendl, R. A. 1999, ApJ, 518, 298 
Putman, M. E., et al. 1998, Nature, 394, 752 

Putman, M. E., Staveley-Smith, L., Freeman, K. C, Gibson, B. 

K., & Barnes, D. G. 2003, ApJ, 586, 170 (P03) 
Recillas-Cruz, E. 1982, MNRAS, 201, 473 

Redman, M. P., Al-Mostafa, Z. A., Meaburn, J., & Bryce, M. 2003, 

MNRAS, 344, 741 
Rich, R. M., Shara, M. M., & Zurek, D. 2001, AJ, 122, 842 
Russel, S. C, &; Dopita, M. A. 1992, ApJ, 384, 508 
Ruzicka, A., Palous, J., & Theis, C. 2007, A&A, 461, 155 
Savitzky, A., & Golay, M. J. E. 1964, Analytical Chemistry 36, 

1627 

Schwarz, U. J., & van Woerden, H. 1974, in Galactic Radio 
Astronomy, eds. F.J. Kerr, & S.C. Simonson III (Dordrecht: 
Reidel), 45 

Sembach, K. R., et al. 2003, ApJS, 146, 165 
Skrutskie, M. F., et al. 2006, AJ, 131, 1163 

Smecker-Hane, T. A., Cole, A. A., Gallagher, J. S., HI, & Stetson, 

P. B. 2002, ApJ, 566, 239 
Smith, A. M., Cornett, R. H., &; Hill, R. S. 1987, ApJ, 320, 609 
Smith, R. C, et al. 1998, PASA, 15, 163 
Snowden, S. L., & Petre, R. 1994, ApJ, 436, L123 
Stanimirovic, S., Dickey, J. M., Krco, M., & Brooks, A. M. 2002, 

ApJ, 576, 773 

Stanimirovic, S., Staveley-Smith, L., Dickey, J. M., Sault, R. J., &; 

Snowden, S. L. 1999, MNRAS, 302, 417 
Staveley-Smith, L., Kim, S., Calabretta, M. R., Haynes, R. F., & 

Kesteven, M. J. 2003, MNRAS, 339, 87 (S03) 
Staveley-Smith, L., Soult, R. J., Hatzidimitriou, D., Kesteven, M. 

J., &; McConnell, D. 1997, MNRAS, 289, 225 
Takakubo, K., &; van Woerden, H. 1966, 

Bull. Astron. Inst. Netherlands, 18, 488 
Townsley, L. K., Broos, P. S., Feigelson, E. D., Garmire, G. P., & 

Getman, K. V. 2006, AJ, 131, 2140 
van den Bergh, S. & Tammann, G. A. 1991, ARA&A, 29, 363 
van der Marel, R. P. 2001, AJ, 122, 1827 

van der Marel, R. P., Alves, D. R., Hardy, E., & Suntzeff, N. B. 

2002, AJ, 124, 2639 (vdM02) 
Wakker, B. P. 2001, ApJS, 136, 463 

Wakker, B. P. 2004, in High- Velocity Clouds, eds. H. van Woerden, 
B. P. Wakker, U. Schwarz, & K. de Boer (Dordrecht: Kluwer), 
25 

Wakker, B. P., Oosterloo, T. A., & Putman, M. E. 2002, AJ, 123, 
1953 

Wannier, P., & Wrixon, G.T. 1972, ApJ, 173, 119 
Wannier, P., Wrixon, G. T., & Wilson, R. W. 1972, A&A, 18, 224 
Weinberg, M. D., & Nikolaev, S. 2001, ApJ, 548, 712 
White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341 
Woosley, S. E., & Weaver, T. A. 1986, ARA&A, 24, 205 
Xu, C, Klein, U., Meinert, D., Wielebinski, R., & Haynes, R. F. 
1992, A&A, 257, 47 



Nidever et al. 



31 



Yamaguchi, R., et al. 2001, PASJ, 53, 985 Zaritsky, D., Harris, J., Thompson, I. B., & Grebel, E. K. 2004, 

Yanny, B., et al. 2003, ApJ, 588, 824 AJ, 128, 1606 

Yoshizawa, A. M., &; Noguchi, M. 2003, MNRAS, 339, 1135 
Young, J. S., & Scoville, N. Z. 1991, ARA&A, 29, 581 
Zaritsky, D., Harris, J., Thompson, I. B., Grebel, E. K., & Massey, 
P. 2002, AJ, 123, 855 



